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ABSTRACT 

The current status is reviewed of theoretical treatments of neutrino-nucleus re- 
actions that are relevant to the detection of astrophysical neutrinos. Various 
nuclear physics aspects involved in the evaluation of the neutrino-nucleus reac- 
tion cross sections are critically surveyed. 
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Chapter 1 Introduction 

Although many topics can come under the title: "Neutrino-Nucleus Reactions", we 
shall be concerned here with neutrino-nucleus reactions that pertain to the terrestrial ob- 
servation of astrophysical neutrinos. The neutrinos play fundamental roles in various astro- 
physical phenomena. The wealth of these phenomena is eruditely described in, e.g., Bahcall 
||^] and Fukugita Q. Obviously, terrestrial experiments to detect these neutrinos are highly 
valuable sources of astrophysical information, and observational neutrino astrophysics is an 
extremely rapidly expanding research field [^. It is useful to distinguish three different 
energy regimes of the astrophysical neutrinos, (i) The low energy region {E,^ ^ 20MeV), 
which includes the solar neutrinos and the lower energy part of the supernova neutrinos, 
(ii) The medium-energy region (20MeV ^ E,, ^ 50 MeV) exemplified by the higher energy 
region of supernova neutrinos, (iii) The high-energy region (50 MeV ^ E^) represented by 
the solar flare neutrinos and atmospheric neutrinos, whose energy can reach as high as 1 ~ 2 
GeV. Neutrino-nucleus reactions offer a great variety of ways to detect these astrophysical 
neutrinos. Generally, the following types of neutrino-nucleus reactions are possible: 

lye + ^A ^ r + ^+^A, (1) 

ue + ^A + ^-^A, (2) 

ue + ^A ^ i^e + ^A, (3) 

i7i + ^A i7(, + ^A. (4) 

In the charged-current (CC) reactions, eqs. the script £ is limited to the electron 

family unless the incident energy is high enough to produce a muon or a tau-lepton, whereas 
the neutral-current (NC) reactions, eqs.(^),(0), can always occur for £ = e, fi and r. The 
final nuclear state A can in general be either bound or unbound (if the latter is allowed 
energetically) , and all these will contribute to inclusive events in which only the final leptons 
are monitored. On the other hand, in a particular type of experiment such as a radiochemical 
experiment, only the contribution of particle-bound final states will be registered. 
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To set the stage, we first give a brief survey of the solar neutrino problem. The Sun 
generates its energy primarily by changing protons into a-particles. The actual reaction 
chain through which p -\- p + p + p — > ^He takes place involves many intermediate nuclear 
reactions, some of which are weak-interaction processes emitting neutrinos. The standard 
solar model (SSM) 0, |] gives a definite prediction of the neutrino flux coming from each of 
these neutrino-emitting reactions. The first observation of the solar neutrinos was achieved 
by Davis and his collaborators through a radiochemical counting of '^^Ar produced in the 
reaction fg + ^"^Cl ^ + '^'^Ar. According to the SSM, the '^'^Cl detector of Davis et aVs 
Homestake facility should register 7.9ib2.6 SNU (1 SNU = 10"^'' events per target atom per 
second), but the observed value is only 2.1ib0.9 SNU This clear deficit of the neutrino 

flux is the solar neutrino problem, which is one of the most fundamental issues in today's 
astrophysics. An independent confirmation of the deficit of the solar neutrino flux has been 



obtained by the Kamiokande experiment 11 , 12 1 , in which the neutrino flux is monitored 



through the leptonic reaction f e + e — > i^g + e occurring in the large water Cerenkov 



counter. The neutrino flux observed by the Kamiokande |12| is 0.46ib0.05(stat) ib0.06(syst) 
times the SSM prediction, corroborating the existence of the solar neutrino problem. One 
thing to be noted here is that the neutrino detection thresholds of the Homestake and 
Kamiokande facilities are rather high. The Kamiokande with -Ethresh = 7.5 MeV is sensitive 
only to the neutrinos coming from the reaction '^B ^Be + e"*" + z^g {Ef^^^ = 15 MeV), 
whereas the chlorine detector is sensitive to the ^B neutrinos and, to a much lesser degree, 
to the neutrinos coming from ^Be + e~ — > ^B + z^g- Since ^B is produced only towards 
the end of the solar thermonuclear reaction chain, the SSM prediction of the neutrino 
flux due to ^B is extremely sensitive to the fine details of solar models. By contrast, the 
neutrino flux due to the primary solar thermonuclear reaction p + p — > -|- e"*" + z^e 
reflects more directly the basic features of the SSM. However, since the maximum neutrino 
energy for this process is small (E^^^ = 0.420 MeV), the detection of the pp- neutrinos 
requires a neutrino absorption reaction with a very low Q-value. The best target so far 



proposed for this purpose is ^^Ga. The threshold of the Uf> + ^^Ga — > e + ^^Ge reaction 
is ^th = 0.2332 MeV, and the SMM prediction for the ^^Ga target is 1321^0 SNU @, of 
which more than half is due to the pp neutrinos. Two groups have recently published the 
results of Ga experiments. The GALLEX collaboration |jl3| reports [83 it 19 (lu) it 8(syst) 
SNU], and the SAGE collaboration ||] gives [58^2! ± 14(syst) SNU], both significantly 
below the SSM value. A wide range of possible solutions to the solar neutrino problem 
have been proposed so far. These "solutions" fall into two categories. In the first category 
(solar-model explanations) a possible deviation from the SSM prediction on the interior 
properties of the Sun is sought. In the second category (particle-physics explanations), the 
observed neutrino flux deficiency is ascribed to unusual properties of the neutrinos. 

Bahcall and Bethe |lf:] have recently carried out a systematic study of solar-model 



explanations. They studied a collection of 1000 precise solar models in which each input 
parameter (the principal nuclear reaction rates, the solar composition, the solar age, and 
the radiative opacity) for each model was drawn randomly from a normal distribution with 
the mean value and standard deviation appropriate to that variable. None of these 1000 
solar models was compatible with the results of the chlorine or the Kamiokande experi- 
ments. Even if the solar models were artificially modified to reproduce a ®B neutrino flux 
in agreement with the Kamiokande experiment, none of these fudged models agreed with 
the chlorine observation. Furthermore, the GALLEX and SAGE experiments were found 
to disagree with the prediction of the solar models by 2o" and 3cr, respectively. Thus, the 
totality of the existing solar neutrino data seems to strongly indicate that the solar neutrino 
problem has at least part of its origin in the hitherto unknown properties of the neutrinos. 

The particle-physics explanations invoke, in one form or another, a mechanism to 
convert electron neutrinos into other particles (neutrino oscillation) so that not all of the 
original u^s produced in the core of the Sun would reach the terrestrial detectors (for a 
review, see e.g. ||l^). The most promising particle-physics solution so far proposed is the 



Mikheyev-Smirnov-Wolfenstein (MSW) effect |18, which can be summarized as follows. 



The three neutrino states, i^e, i'^ and are assumed to have different masses, with the 
electron neutrino the Hghtest. In the solar core with high electron densities, the neutrinos 
acquire effective masses through the neutrino-electron interactions. The neutral-current 
interaction simply gives a universal shift of the neutrino masses, but the charged-current 
interaction, which in the absence of muons and r-leptons is operative only for i/g, increases 
the rriu^ selectively so that, for sufficiently high electron densities, the effective electron 
neutrino mass can become heavier than those of the other neutrinos. If the Hamiltonian 
involves a piece that causes lepton-flavor mixing, then in a critical region where m'^ ~ m'^, 
or m^^ ~ m'^, it is possible for Ug to transform into the other neutrino species via reso- 



nant transitions. This phenomenon is called the MSW effect ||l^, 19]. The parameters 
characterizing the MSW mechanism are the neutrino mass differences in free space and the 
strengths of flavor mixings. By adjusting these parameters one can vary the region of the 
electron neutrino spectrum that is changed into the other neutrinos through the resonant 
transitions. The contemporary literature abounds in articles that deal with the determina- 
tion of the parameters characterizing the MSW model or its variants. At present, the bare 
neutrino mass difference of about 10^^ to 10^^ eV with appropriate corresponding mixing 
strengths is compatible with the existing data. However, this does not constitute unam- 
biguous evidence that the MSW mechanism is operative in the Sun; the direct confirmation 
of the MSW effect needs an additional set of experiments. Neutrino-nucleus reactions are 
of great importance in this context, since in some favorable cases they allow us to identify 
various types of neutrino reactions individually. Of particular significance is the detection 
of neutral-current (NC) reactions, eq.(^); since the NC reaction rate is the same regardless 
of neutrino species, a measurement of the NC reactions gives the bolometric flux of all neu- 
trino species |20, 21, 22]. Therefore, by combining information on the NC reactions with 



that on the CC reactions, one can obtain definitive information on neutrino oscillation, or 
more specifically, on the MSW effect. In the planned Sudbury Neutrino Observatory (SNO) 



]23, 24], which uses a 1,000-tonne, heavy-water Cerenkov counter, one aims at the separate 



registration of the CC and NC reactions on the deuteron with the detection threshold as low 
as Ely =5 MeV. The data obtained at the SNO will provide unambiguous information on 
the neutrino flux for individual flavors and consequently on the neutrino oscillation. There 
are other nuclear targets which are important in this context, and we will discuss them later 
in the text. 

The neutrino-nucleus reactions are not the only place where neutral-current interac- 
tions appear. The leptonic reactions: 

i^£ + e~ i^e + e~ , (5) 
1^1 + e~ ^ Ui + e' , (6) 

where i = e, or r, do involve the neutral-current contributions. More specifically, for £ = e, 
the reactions can occur through both the charged and neutral currents, whereas for i = fj, 
and T only the neutral current will participate. Therefore, water Cerenkov counters and 
liquid scintillation counters in which these leptonic reactions provide main signals do count 
both neutral-current (NC) and charge-current (CC) events. It is to be noted, however, that 
the NC and CC events cannot be registered separately by monitoring the leptonic reactions, 
eqs.(^),(^). This situation makes the use of nuclear targets particularly attractive. 

We now turn to astrophysical neutrinos other than the solar neutrinos. The large 
water Cerenkov counters at the Kamiokande and 1MB allowed the first observation of the 
supernova neutrinos coming from SN1987A. For this spectacular success it was crucial that 
these facilities had sufficiently high time resolutions to clearly define supernova-related neu- 
trino events. Since all the existing nuclear-target detectors use the radio-chemical counting 
method, their time resolution is not sufficient to allow the monitoring of supernova neutri- 
nos. In the near future, however, counter experiments involving nuclear targets will become 
possible. The above-mentioned Sudbury Neutrino Observatory (SNO) [p^, 24 1 is a notable 



example. As will be discussed in more detail later, the SNO's ability to detect the NC 
reactions with high time-resolution can be important in the context of the supernova neu- 



trinos as well. We will also discuss some other nuclear targets which can be very useful in 
connection with supernova neutrino detection. 

In considering the high-energy astrophysical neutrinos [group (iii) in the above- 
described classification], it is important to realize that, as the incident neutrino energy 
increases, the nuclear reactions become progressively more important as compared to the 
leptonic processes. This implies that even the water Cerenkov counter and liquid scin- 
tillation counter, which for low-energy neutrinos primarily detect the electrons from the 
leptonic reactions, will start detecting the nuclear reactions |^]. As will be discussed in 
the following, this offers yet another interesting use of these facilities. A notable example 
is the water Cerenkov counter measurement of the ratio of the muonic neutrino flux to the 
electronic neutrino flux in the atmospheric neutrinos. 

One of the fundamental prerequisites for nuclear targets envisaged as neutrino de- 
tectors is that reliable estimates of the relevant reaction cross sections be available. This 
presents an interesting and challenging task to nuclear physicists. f\ The great variety of 
relevant nuclear targets and the wide energy range of the astrophysical neutrinos to be con- 
sidered necessitate the use of many different approaches that can be found in the nuclear 
physics arsenal. In this review, we survey the current status of theoretical attempts at 
obtaining reliable estimates of the neutrino- nucleus reactions of astrophysical interest. To 
highlight the basic aspects of various approaches so far used and their possible limitations, 
we shall organize our topics according to the calculational methods involved. We describe 
in chapter 2 the direct microscopic calculation which is relevant to the lightest nuclei. In 

chapter 3 we discuss a model-independent determination of the cross sections based on 
^The cross sections of the leptonic processes, eqs.(|^,(^), by contrast are accurately known from the 
standard electroweak theory. There is in fact one particular type of neutrino-nucleus reaction for which the 
cross section can be determined mode-independently; for the coherent scattering v + yl(gnd) v + A(gnd), 
the cross section is determined by the number of nucleons in A. The coherent process plays an important 
role in dark matter search as well as in certain types of solar neutrino experiments ^ , but we 

will not discuss the coherent scattering here. For a review of this topic, see M . 
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the elementary-particle treatment. Chapter 4 deals with the empirical effective operator 
method. The use of (p,n) reactions as calibrators of the Gamow- Teller matrix elements is 
surveyed in chapter 5. Chapters 6 and 7 are concerned with high-energy neutrino-nucleus 
reactions. The Fermi gas model and its possible refinement are discussed in chapter 6, while 
the treatment of semi-inclusive reaction is addressed in chapter 7. Additional remarks are 
made in chapter 8. 
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Chapter 2 



"Ab initio" calculations based on realistic wave functions 



In the so-called realistic description of nuclei one works with the nuclear Hamiltonian 



where A is the mass number, Tj is the nucleon kinetic energy, and Vij is the realistic nucleon- 
nucleon potential. Although Hf^f is called the realistic Hamiltonian due to the fact that Vij 
reproduces the observed properties of the two-nucleon systems reasonably well, is an 
effective Hamiltonian of highly operational nature. To arrive at starting from the fun- 
damental QCD picture is far from trivial. First, the dynamical degrees of freedom of quarks 
and gluons need to be translated into the effective degrees of freedom of hadrons. Second, 
the Hilbert space of hadrons must be truncated down to those of non-relativistic nucleons 
interacting via potentials. Despite these fundamental problems, it is generally believed that 
solutions for the Schrodinger equation with the Hamiltonian give reasonably reliable 
nuclear wave functions, and it is usually considered as an ab initio approach to calculate 
matrix elements of nuclear electromagnetic or weak-interaction processes with the use of 
these wave functions along with the standard non-relativistic reduction of the single-nucleon 
responses to the relevant currents Two major problems encountered in this ab ini- 

tio approach are as follows. First, it is in general extremely difficult to obtain an exact 
solution of the ^-body Schrodinger equation Hn'^{1, 2, ... A) = 2, . . . A), and hence 

one is usually forced to work with rather drastically truncated model wave functions ^'o; 
shell-model wave functions are the most commonly used and most successful example. Now, 
if the matrix element of a nucleon observable On is calculated using model wave functions 
and ^l, one expects < ^'•^ | Oat | >/< ^1\On\% >■ This difference is called core- 
polarization effect. Another point is that eliminating all hadronic degrees of freedom but 
those of the nucleons is an approximation and, even if one can calculate < ^'■^IC'tvI^* > 
exactly, it may still differ from the corresponding true matrix element. One usually ac- 
counts for this difference by adding to On many-body operators O' , which represent the 



A 



A 




(7) 
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effects of the eliminated hadronic degrees of freedom. These additional terms are called the 
exchange-current operators. 

Now, to calculate the core-polarization and exchange-current effects in complex nu- 
clei is one of the most challenging problems in nuclear physics and, although many impor- 
tant achievements have so far been made in this domain, the problem is far from solved. 
However, nuclei with small mass numbers {A < 4) offer important exceptions. For these 
lightest nuclei one can obtain sufficiently "realistic" wave functions either by directly solv- 
ing the Schrodinger equation or through variational calculations, eliminating thereby the 
core-polarization problem. For these systems, therefore, one can carry out ab initio calcu- 
lations of the transition matrix elements of electromagnetic and weak-interaction processes, 
provided a reasonably reliable method to incorporate the exchange-current effects exists. 

2.1. u — d reactions 

The A=2 systems are the simplest nuclei, for which one can obtain reasonably reliable 
wave functions corresponding to the realistic Hamiltonian eq.(0). This allows us to esti- 
mate with high accuracy the cross sections of the neutrino-deuteron reactions. These cross 
sections are important for deriving useful astrophysical information from the experimental 



data that will soon become available at the Sudbury Neutrino Observatory (SNO) 24 1. 
The neutrino-deuteron reactions that are relevant to the SNO are: 

v + d^u' + n+p, (8) 

u + d^v + n + p. (9) 

Ve + d e-+p + p, (10) 

r'e + d — > + n + n, (11) 

As mentioned earlier, one of the great advantages of the SNO is its capability to detect 

both the neutral- and charged-current reactions and record them separately, and this is 



crucially important [22| for proving that the MSW effect |18| is indeed a solution for the 
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solar neutrino problem ||T|. 



The earlier calculations 33, |3J] were mainly meant for low-energy solar neutrinos 
and therefore included only s-wave for the final N-N relative motion. Extended calculations 
that take account of higher partial waves and, correspondingly, the contributions of the 
forbidden-type transitions were carried out by Ying, Haxton and Henley [^5| (to be referred 
to as YHHl) and by Tatara, Kohyama and Kubodera (to be referred to as TKK) up 
to El, = 170 MeV. Unfortunately, the cross sections obtained in these two detailed calcu- 
lations show alarmingly large discrepancies for higher energies; the YHHl results become 
significantly larger than those of TKK around E^, ~ 20 MeV and the cross section ratios 
reach ~ 2 towards Ei, = 170 MeV. To clarify the origin(s) of these discrepancies, Doi and 
Kubodera ||3^ (to be referred to as DKl) carried out a systematic comparison of the for- 
malisms used in YHHl and TKK. The numerical results of the independent calculation in 



DKl |37| suggested possible numerical errors in YHHl, and this suggestion was confirmed 
by the authors of YHHl. The results of a revised estimation by Ying, Haxton and Henley 
|l38|| (to be referred to as YHH2) are in essential agreement with those of TKK and DKl. 

Beyond the ~20 % level precision, however, there still exist discrepancies among the 
results of TKK, DKl and YHH2. The cross sections in TKK begin to be systematically 
smaller than the corresponding numbers in DKl and YHH2 as Ey surpasses ~ 100 MeV, 
and the differences for some channels reach ~20 % towards Ei, = 170 MeV. It has been 
demonstrated in DKl that this discrepancy arises because TKK dropped the term quadratic 
in , the space component of the vector current. At low energies, the contribution of 
is known to be much less important than that of the time-component pv- TKK, assuming 
that this behavior would persist for higher energies, ignored the quadratic term | that 
appears in squaring the transition matrix element to obtain the cross section. According 
to DKl, however, the contribution of is comparable to or even larger than that of pv 
for E,y ^ 10 MeV, and ignoring | leads to an underestimation of the cross sections by 
up to ~20 %. Meanwhile, cr{i'e -|- d — > e+ + n + n) of DKl becomes progressively larger 
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than the corresponding results in TKK as Ei, decreases from 7 MeV; e.g., at = 6 MeV, 
a{i'e + d — > 6+ + n + n)DKi/c(^'e + d — > e"*" + n + n)TKK = 1-07. (At this low energy 
the contribution of | ignored in TKK is completely negligible.) A recent study by Doi 
and Kubodera (to be referred to as DK2) indicates that this behavior is caused by the 
approximate Q- value expression used for this channel in DKl. As far as comparison with 
YHH2 is concerned, we remark that the exchange-currents effects are altogether ignored in 



YHH2. This is estimated [36| to lead to a general underestimation of the cross sections by 



up to ~5 %. Thus each of the published works, TKK [||], DKl and YHH2 |8|, has 
some point (s) to be improved upon. 



Kohyama and Kubodera |40| have extended the calculation of Tatara et al. |3£] by 
including the contributions of the terms. The results of Kohyama and Kubodera 



1 40 1 and those of DK2 [39| have been found to agree with each other within 2 %. Since 
these two calculations are independent, both in the calculational methods and in choosing 
various input parameters, it is reassuring that their results agree to this degree. (According 
to TKK, different choices of the nuclear potential from among realistic nucleon-nucleon 
interactions affects the cross sections by up to 1 %.) We tabulate in table 1 in Appendix 
the cross sections for the reactions eqs.(8)-(ll) obtained by Kohyama and Kubodera pO|| ; 
5 % uncertainties are attached to the cross sections. Table 1 represents the best available 
estimates of the neutrino-deuterium cross sections up to Ei, = 170 MeV. 

One generally expects that the ratios Ru = (T(t'e + c? +p+p) / (T[i'+d — > v' +n+p) 
and Ry = (T{D(, + d — > e'^ + n + n)/a{h' + d — > v + n+p) are much less sensitive to the nuclear 
models than the absolute cross sections themselves. Indeed, the R^, and Rp obtained by 



Kohyama and Kubodera |4C] agree with the corresponding values given in DK2 |39] within 
0.2 %. Therefore, it seems reasonably safe to assign 1 % uncertainties to the calculated 
values of Ry and 

The average cross section for the electronic neutrinos from the /i"*" decay calculated 



using the results of |40| is: < a{i'e + d— >e + p + p) >/i+-decay= 5.4 x 10 cm^. This 



IS 
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consistent with the observed value a = (5.2 it 1.8) x 10 cm^ pl[| . 



Apart from testing the MSW effect for the solar neutrinos, Bahcall et al. [34| con- 
sidered the merit of the SNO in determining the masses of heavy-flavor neutrinos through 
arrival-time-delay measurements of supernova neutrinos. Some other possible uses of the 



SNO for solar and supernova neutrinos were discussed by Balantekin and Loreti [42|. The 
SNO can also be used to study neutrinos of extremely high energies such as solar-flare 
and atmospheric neutrinos. The above described "a6 initio" calculation probably can be 
extended up to E^, ~ several hundred MeV, but its reliability will diminish rather signifi- 
cantly for these high energies. At yet higher energies, the results are expected to approach 
those of the high-energy approximation to be discussed in chapter 6. 

2.2. Neutrino-production reaction on the lightest nuclei 

Although our review is primarily concerned with neutrino- nucleus reactions, this is 
a good place to discuss the latest development in the theoretical treatment of neutrino- 
production reactions on the lightest nuclei. As mentioned earlier, the pp reaction 

p + p^d + e+ + Ue, (12) 

is the primary solar thermonuclear reaction and produces a dominant part of the solar 
neutrino flux Q. The "/lep" reaction 

^He + p ^ ^He + e+ + i/g , (13) 

contributes only a tiny portion of the solar neutrino flux but the maximum neutrino 
energy of the hep process {E^^^ = 18.77 MeV) is the largest of all the solar neutrinos. 



Carlson et al. |43| have carried out a detailed ab initio calculation of the cross sections 
for these reactions, using as a nucleon-nucleon interaction the Argonne V14 potential. To 
reduce the uncertainty in the axial exchange current operator, its matrix element has been 
adjusted so as to reproduce the measured Gamow- Teller matrix element for tritium /3-decay. 
For the pp reaction, the calculated astrophysical S factor and its derivative at zero energy 
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are: S{E = 0) = 4.00 x lO'^^ MeV b, and dS{E)/dE\E=o = 4.07(1 ± 0.051) x lO^^s MeV 
b; these are close to the values quoted in Carlson et al. find that exchange-current 
contribution enhances the capture rate by 1.5 %. This is considerably smaller than the 



enhancement found in earlier works |44, 45|. In the absence of Coulomb distortion the 
transition matrix element of the pp reaction is just the hermitian conjugate of that of 
the fg + d— > e~ + p + p reaction. For the thermal pp reaction, however, the Coulomb 
repulsion between the two protons is extremely important and responsible for the reduction 
of exchange-current contribution as compared with the typical values found for p^, + d ^ 
e~ + P + p reaction fsH] . 



As for the ^He(p, e^i^)'*He reaction rate, Carlson et al. |43] obtained S{E = 0) = 
1.3 X 10^^^ MeV b, which is smaller by a factor of ~ 6 than the values obtained earlier by 



Werntz and Brennan [[iq ] and by Tegner and Bargholtz [47|. This change is caused primarily 
by the exchange current effect, which cancels the impulse approximation contribution. This 
is one of the most dramatic examples of the exchange-current effect. 
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Chapter 3 Elementary-particle treatment (EPT) 

The elementary-particle treatment (EPT) was first introduced by Kim and Primakoff 



|48| and by Fujii and Yamaguchi [49|. In EPT, instead of describing nuclei in terms of nu- 
cleons or other constituents, one treats nuclei as "elementary" particles with given quantum 
numbers. A transition matrix element for a given process is parametrized in terms of the 
nuclear form factors solely based on the transformation properties of the relevant current 
and nuclear states |5^. Insofar as all of these nuclear form factors can be deduced 

empirically from experimental data, one can make a totally mo del- independent prediction 
on every observable for the given transition. In fact, this requirement is an extremely strin- 
gent one, and so far the EPT calculation has been carried out completely only for the triad 

of -12c -12N. 

This particular A=12 case, however, deserves special attention. One reason is that 
the 12c nucleus, which is abundantly contained in ordinary liquid scintillators, can be a very 
important nuclear target in neutrino astrophysics. Although the threshold energies of the 
neutrino reactions on 12 C are too high for solar neutrino detection, they are low enough for 
observing energetic neutrinos originating from stellar collapses. The neutrino energies up to 
El, ~ 100 MeV are expected to be of relevance for stellar collapse neutrinos. Another reason 
is provided by the intensive programs of beam-dump neutrino experiments at Los Alamos 
|52, 53, 55| and the Rutherford Laboratory 57|. Some of these experiments directly 



aim at observing the neutrino -12 C reactions. In other particle-physics experiments 
main interest is in the — e scattering but here also, to eliminate the dominant background 
due to the neutrino -12C reaction, one needs to know its cross sections with reasonable 
accuracy. 

The processes of relevance here are the superallowed neutral-current (NC) reactions 

i2C(0+;gnd) ^ i2c*(l+;15.1MeV), (14) 

i2C(0+;gnd) ^ P'+ i2c*(l+;15.1MeV), (15) 
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and the superallowed charged-current (CC) reactions 

^2c(o+;gnd) ^ e" + ^2N(i+;gnd), (16) 
P+ i2C(0+;gnd) ^ e+ + ^2B(;L+.gnd-)^ (;L7) 

The final nuclear states of the above reactions form a triad of = 1+,T = 1 states. The 
existence of these superallowed transitions is favorable for simultaneous monitoring of the 
NC- and CC-reactions. (The elastic scattering leading to ^^C(0^;gnd) is unobservable in 
normal circumstances.) Fukugita, Kohyama and Kubodera [^ ] showed that EPT allows the 
completely model-independent determination of the cross sections for the reactions eqs.(p^- 
(|l3) for ^ 100 MeV. A similar treatment was presented by Mintz and Pourkaviani 
p9|, pO|]. Although the use of EPT is at present limited to the transitions listed above, this 



is not an obstacle for monitoring these exclusive processes, as discussed in |58]. Here we 
summarize the calculation of Fukugita et al. |^] . 

The effective Hamiltonian responsible for the NC reactions is given by 

H,jf = GF/V2JlL% (18) 

where Gp is the Fermi constant, Lj^ = 'il>u^^{l — j5)ip,y, and = (1 — 2sin^ 0w)y^ + 
with 9w being the Weinberg angle. (The isoscalar neutral current, which cannot cause 
nuclear excitations, has been dropped.) The effective Hamiltonian for the CC reactions is 
given by 

Heff = [Gf cos 9JV2] jlt^L^^\ (19) 

where 9c is the Cabibbo angle; = '4>elii{'^-l^)tpv, and L^^^ = L^T^^ J^^^ = V^i^^+A^if^ . 
The most general form of the matrix elements of the vector and axial-vector currents for 
the NC reactions is 



C(l+);p,|4(0)|i2c(0+);^,, >= ^[^^FAq^) + q^{C.q)^^ +Q,{C-q)^ 
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where q = Pb — Pa, Q = Pb + Pa, ^ is the polarization vector of the spin-1 nucleus, and 
M = (Ma + Mf,)/2 is the average of the initial and final nuclear masses. All information 
on nuclear dynamics is contained in the nuclear form factors, FM{(f'), FA{q^), Fp{q'^), and 
Fxiq'^). Here the nuclear form factors are classified in the "cartesian" representation. An 



alternative formalism is described in |5C, 62]. For the CC reactions the Wigner-Eckart 
theorem in isospin space leads to 

N(l+);p,|4+)(0)|i2c(0+);p, >= C,FA{q') + q,{C ■ q)^^ + Q,{i ■ Q)^^- (23) 

Now, the contribution of the Fp term to the observables is proportional to (lepton mass)^ 
and therefore can be dropped in the present energy regime {E,, ^ 100 MeV) for which muon- 
production can be ignored. For the remaining nuclear form factors one can first consider 
their values at q'^ = 0, and examine their g^-dependence later. Define Fm = FM{q^ = 0), 
Fa = FA{q^ = 0), and Ft = Frp{q'^ = 0). Fm can be determined from the width of the 
7-decay ^'^C*{1+; 15.1MeV) ^ 12^(0+; gnd), and the CVC. From F^^p = 38.5 ± 0.8 MeV 
p3| , we obtain Fm = (1.516 ±0.016) x 10^'^ MeV~^. Fa and Ft can be determined from the 
/?+ decay i2N(l+;gnd) ^ ^^Q^Q+jgnd) and the /?" decay i2B(l+;gnd) ^ i2Q(o+;gnd). 
The /i- value of the /3-decay is written as 

ft = ''-f:^ (24) 

■' ml{GFCOse,YFl{l + lEoFT/FAy ^ ' 

with Eq the maximum /3-particle energy. The measurement of the /?-ray angular distribution 
from aligned ^^B and ^^N parent nuclei gives the angular correlation quantities a±, which 



are related to Ft as 65, 36 1. 



a+ + a- = -^FT/FA. (25) 

The experimental values |6^, |8|, a_ = (0.02 ± 0.12) x lO-^MeV"^ and a+ = -(2.72 ± 
0.27) X 10"3MeV"\ give Ft/ Fa = (2.01 ± 0.23) x lO^^MeV"^ Then the contribution 
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of the Ft term in eq.(|2|) is {2/3)Eo{Ft/Fa) = 0.02 for Eq = 15 MeV. The importance 
of this term relative to 1 being smaller than typical isospin-symmetry breaking effects 
within nuclei {aZ ~ 6%), the Ft term in eq.(]2^) can be ignored. Experimentally | |63| , 
log(/i) = 4.067 ± 0.002 for ^ ^^C, and log(/t) = 4.120 ± 0.003 for ^ ^^C. The 
difference between {ft)^ and (ft)^ gives a measure of the isospin symmetry breaking quoted 
above (insofar as there is no second-class currents). Using ft = 12444 it 830 s, which is 
the average of (ft)^ and (ft)^ with the error corresponding to their difference, we obtain 
Fa = 0.711 ±0.024. 

As E,^ grows, one must consider the g^-dependence of the nuclear form factors. The 
total cross section for the NC reactions including the g^-dependence is given by 

a = aoF^E^d, (26) 
(To = 3G^/2^ = 2.53 X 10-^^cmV(MeV)2, (27) 
I = i|'^dz/(q2)[(l_£)±l(i?, + i^^,)(l_2sin20iyTy)(l-z)^ 
+ 1[E,>E,{1 - z") + (1 - z)q\l - 2sin2 Owfi^f 
_|AM(l+„g+l,^(g^„. ,28) 

where the upper (lower) sign refers to the neutrino (anti-neutrino) reaction; z = cos 9 
with 9 being the neutrino scattering angle; E^' = E^ — AM with AM = 15.1 MeV; 

= [^Aiq^)/ Fa\^ with = (AM)^ — q^. The cross sections for the CC reactions are 
given by 

a = 2(70 cos 9cPeEeF^{Zf; E^)fIi, (29) 

where the upper (lower) sign refers to the neutrino (anti-neutrino) reaction. Here / is as 
defined previously except the replacements E^i — > i?e, and (1 — sin^ ^vk) 1; F^{Zf\Ee) 
denotes the Fermi function for the Coulomb correction. Ef> = E^ — AM, with AM = 16.827 
MeV (13.880 MeV) for {^^C ^^B). The above expression has been obtained 

by assuming that the q^-dependences of FM{q^) and FT{q^) are the same as that of FA{q'^)- 
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From an EPT analysis of the "^^C + /u~ ^ ^^B + process, Nozawa et al. showed that 
Fa{q^) / Fa = FM{q^) / Fm holds within a 10 % accuracy at least up to = 0.74m^. For 
Fxiq'^) no direct experimental information is available but, since its contribution is minor, 
its precise form has little significance. Now, /(q^) — [FM{q'^) / Fm^ can be determined 
model-independently from FM{q'^) obtained from inelastic electron scattering: e + ^^C — > 
e'+ ^2c*(l+;15.1MeV). According to p|. 



fiq') = [1 - (1 - p)mq\)']' exp[-i(6|q|)2], (30) 

where b = 1.881 ± 0.053 fm, and p = 0.23 (1 ± 0.1). 

Using the EPT summarized above, Fukugita et al. [58| were able to determine 



the cross sections for eqs. (p^)-(|T7|) up to E^, ~ 100 MeV within 12% accuracies; 
these uncertainties reflect the ambiguities in the input data and isospin symmetry break- 



ing in nuclei. Similar results were obtained by Mintz and Pourkaviani S0|. As in 
the A=2 case discussed earlier, the ratio of the NC cross section to the CC cross sec- 
tion is expected to be less affected by the nuclear physics ambiguities. When we vary 
the nuclear form factors within ~10 % (which are typical errors given in ref. |^8|), R = 
a[^^C{ue, u'^y^C*{l+; 15.1MeY)]/a[^^C{iye, e~y^N{gnd)] is found to change by less than ~2 
%. 

We now remark on another semi-empirical method developed by Walecka and Don- 
nelly to estimate the matrix elements of weak-interaction transitions in light nuclei. In 



the Walecka-Donnelly (W-D) method |71, one assumes that the weak-interaction 

transition operator can be approximated by a sum of single-particle multipole operators 
that arise in the impulse approximation and that the electomagnetic data can be used to 
place constraints on the nuclear matrix elements of these operators. Because of these basic 
assumptions it should be obvious that, contrary to the often-made claim, the W-D method 
is not a model-independent approach. Apart from the fundamental assumption pertaining 
to the impulse approximation, an important question is to what extent one can identify the 
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effective matrix element of a single- nucleoli operator On appearing in the weak-interaction 
process with that of the On that appears in the electromagnetic process. Chemtob and 



Rho [75 1 emphasized that the exchange-currents for the electromagnetic interaction can be 



drastically different from those of the weak processes. This point has been reemphasized 



in |43|. The EPT is free from these problems and hence, from a formal point of view, 
definitely a preferable approach. Having said this, we must quickly add that EPT in its 
rigorous form can be used at present only for the very limited case of A=12. By contrast, 
the W-D method, owing to the additional assumptions made on the possible form of the 
transition operators, can be applied to many cases where the lack of necessary experimental 
information renders EPT at moment powerless. Thus, in practice, the W-D method is a 
very useful approach insofar as some independent information is available to check its relia- 
bility for given individual cases. The use of experimental data to constrain the parameters 
involved in the W-D method is in general expected to reduce nuclear-model dependence 
significantly, but the remaining possible uncertainties should be kept in mind. 



Donnelly and Peccei |21] used the W-D method to estimate the cross sections for the 
NC reactions eqs. (|T^) , (p^) . While a good agreement is seen between Donnelly and Peccei's 
results and those of Fukugita et al. js^ at lower incident energies, disagreement becomes 
apparent for ^ 25 MeV, and it amounts to 60 % at ~ 50 MeV. Furthermore, while 
Fukugita et al. report an appreciable difference between ^'^C{u,u'y^C* and ^^0(2^, z^')^^C* 
at higher energies due to the interference term FaFm, this behavior is not reproduced in 



|21|. These discrepancies seem to signal the breakdown of some assumptions involved in the 



W-D method p^] as applied to higher incident neutrino energies. On the other hand, for low 



neutrino energies, the EPT results justify the W-D method |71], in which the leading-order 
impulse-approximation transition matrix element < tct > in the weak interaction process is 
deduced from the observed Ml strength. For ^^C(z^e) e~)^^N the estimate of Fukugita et al. 



| 58| is significantly smaller than the published result of the W-D method |7^; specifically. 



for the cross section averaged over the Ve spectrum due to the /i-decay, the EPT gives 
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< cr[^2C(i/e,e~)^^N(gnd)] >= 9.2 x lO^^^cm^, which is about 25 % lower than that of 



Donnelly |73|. However, a more recent result of the W-D method agrees with the EPT 
result. 

The KA rlsruhe-Rutherford Medium Energy Neutrino (KARMEN) experiment has 



been intensively pursued at the spallation facility of the Rutherford Laboratory |£6|, |76|. 
The pulsed beam dump neutrino source provides monoenergetic i/^ with Eu = 29.8 MeV 
from vr^ decay at rest as well as z^e and with energies up to 52.8 MeV coming from 



the subsequent fi^ decay. A 56 t liquid scintillation calorimeter [76| is used as massive 



live target of ^^C. The NC reaction eq.(^) has been observed |77| by monitoring the 7 
emission back to the ground state {E^ = 15.1 MeV, B.R = 94 %). Experimentally |0, < 
o-[i2C(i/e, Ky^C*{l+; 15.1 MeV)] >''^p= [9.5±1.8(stat)±1.35(syst)] x lO'^^ ^^2^ ^j^^^g ^ ^> 
represents the average over the Ve spectrum coming from the decay. This agrees well with 
theoretical estimates < a[^'^C{ue,u'^)^^C*{l+;15.1MeY)] >thcor^ (g g _ ;^q 3) ^ iQ-^a cm^ 
|58, |6^, The ^^C(z^e; e~)^^N(gnd) reaction can be identified by taking the delayed 

coincidence of the recoil e" and the e"*" from the ^^N /3-decay. The KARMEN measurement 
M, 0] gives < a[^^C{i'e,e-)^^l<l{gnd)] >'^^p= [8.0 ± 0.75(stat) ± 0.75(syst)] x lO'^'^cm^ 



which is consistent with the earlier Los Alamos result ^3|: <a[^'^C{i'e, e )^^N(gnd)] >''^p= 
[10.5 lb l.O(stat) lb l.O(syst)] x 10^^^ cm^. These experimental values are in good agreement 
with recent theoretical estimates <a[^^C{i'e,e~y^Nignd)] >thcor^ Q _ g_4^) ^ iq-^s 



given in |5|, |5|] . 

EPT calculations on the u-^'^C reactions for Ei, ^ 100 MeV were carried out by 
several authors. In particular, Mintz and Pourkaviani [^] extended the EPT calculation 
up to = 3 GeV. The reliability of such extension, however, is not well established 
since the calculation involves the nuclear form factors at very high q'^, which are not known 
empirically except Fm{q'^)- There have also been attempts to apply EPT to nuclear systems 



other than the ^^C target or to processes [59|. Here again, the problem is that 



there is not enough data to determine all relevant nuclear form factors empirically. The 
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lack of experimental information necessitates the introduction of a number of additional 
theoretical assumptions, which are usually motivated by the impulse approximation results. 
This ad hoc procedure largely nullifies the predictive power of the original EPT. 

Koetke et al. [E4|, using the in-flight pion decay neutrino source at LAMPF, measured 



the cross sections for "^"^(^{u^, fi^)X and ^^C(f^, /i^)"'^^N(gnd) at an average neutrino energy 
of 202 MeV. The cross section for the latter exclusive process is (TpC(i/^, /z~)"'^^N(gnd)]''^P = 
(1.7 ± 0.8 ± 0.3) X lO^'^^cm^. This is considerably higher than 0.8 x 10""^° cm^ obtained 
in the W-D method |8^, |2[ and an "EPT" estimate of 0.7 x 10"^° cm^ |9[. Here we 
have used the quotation mark since, as discussed above, a completely model-independent 
EPT analysis is not feasible in this case. According to Mintz and Pourkaviani one 
can increase the predicted value sufficiently by enhancing the pseudoscalar form factor fp 
Although it is generally believed that the exchange-current effect quenches fp rather than 
enhances, there may be other competing nuclear effects. (For recent attempts to determine 



the effective fp inside nuclei using the radiative ^-capture, see ref.[33| for experiment and 
ref.|84] for theory.) The inclusive reaction ^'^C{u^j_, n^)X will be discussed later. 
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Chapter 4 Empirical effective operator method (EEOM) 

Going back to the "realistic" nuclear Hamiltonian Hj^ eq. (|^, we recall that, even 
if we forgo the basic problem of deriving from the fundamental QCD and decide to 
be content with the operational usefulness of Hn, it is in general prohibitively difficult to 
determine true eigenstates ^''s of and calculate the nuclear matrix elements of an ob- 
servable O using ^''s. This is true even when we limit ourselves to the nucleon-only regime, 
ignoring the exchange-current problem. One therefore usually splits the space subtended 
by ^'s into two parts, a relatively manageable model space P and the remainder Q, and 
calculate observables using wave functions ^'o's that belong to P. The phenomenologi- 
cal success of the shell model suggests choosing as P the space of the lowest shell-model 
configurations; if necessary, slightly excited configurations can be included in P. Once 
the space of 'f's is truncated to that of ^'o's, Hn and O must undergo corresponding 
transformations, Hn Hn and O ^ O so that < "^'IHn]^ >=< ^oI-^a^I^o > and 
< ^''|C'|^ >=< ^'qIC'I'I'o > should hold. The nucleon-nucleon interactions that feature 
in Hn are called the effective interactions, a somewhat historical terminology since from 
today's viewpoint the original Hf^f itself is a highly "effective" entity. Similarly, O is called 
the effective operator. The difference between O and O represents the core-polarization 
effect. The formal framework to obtain the effective operators by incorporating effects of 



the eliminated space Q does exist (for review, see [85, q6|), but applying this formalism to 
actual complex nuclei encounters many practical difficulties. Exceptionally, in the simplest 
cases of closed-shell nuclei plus or minus one nucleon, explicit calculations can be carried 
out to lowest orders of perturbation, as exemplified by the extensive work of Tower and 



Khanna |92, 93 1 and the University of Tokyo group (Arima, Shimizu, Hyuga et at) 
These authors calculated the magnetic moments. Ml transition strengths -B(Ml), and the 
Gamow- Teller strengths -B(GT), taking into account up to second-order core polarization 
effects. Furthermore, the exchange-currents effects were carefully estimated and included 
in the final expressions of their effective operators. The results of these highly elaborate 
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calculations are thoroughly documented in p^ ]. 

Parallel to these theoretical attempts to derive the effective operators for the special 
cases, a useful semi-empirical approach which can cover a wider variety of nuclei has been 
developed by Wilkinson, Brown, Warburton and Wildenthal [^, ||, |9^, |9l| . We shall 
refer to this approach as the empirical effective operator method (EEOM). As will be 
explained below, EEOM finds a useful application in estimating the cross sections of low- 
energy neutrino-nucleus reactions for nuclei up to the sd-shell. We outline EEOM, taking 



the Gamow- Teller (GT) operator as an example |91|. 

Consider a GT transition \i >— > \f > and suppose we are given reasonably realistic 
nuclear wavefunctions for \i > and \ f >■ By "realistic" we usually mean the best available 
shell-model wavefunctions, which, as today's standard demands, are eigenvectors of a full 
matrix of Hj^ covering all possible configurations within a given major shell. The shell- 
model GT matrix element is given by 

A 

<GT> = <f;Jf,Tf\\\Y,gATCT{k)\\\i;J,,T,>, 

k=l 

= VeY^^fj <f\\9A(r\\j>, (31) 

where the summation over j and j' goes over all the single-particle orbits contained in 
the model space P, and the one-body-transition density Dj/j is obtained from shell-model 
wavefunctions for an operator with rank in ordinary (isospin) space of A J (AT) via 

n < J/,r;|||[at(/)xa(i)]-^'^-^^|||J.,r.> 
[(2AJ + l)(2Ar + l)]V2 • 

As stated, even if the model space P is reasonably realistic, there are two significant sources 
of corrections to < GT >, one from core polarization and the other from the exchange- 
current effect. If P is large enough, the core polarization effect arises from rather highly 
excited states, and consequently its dependence on |i > and \f > will be weak. The same 
should hold for the exchange-current effect, which comes from the eliminated hadronic de- 
grees freedom. In EEMO, therefore, one simulates these effects by replacing < /Hfi'yicrllj > 
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in eq. ( pjl| ) with an effective matrix element 

<j'\\9fcT + gfAl + gfA[Y2 X ^](^)|U >, (33) 

while keeping Djij unchanged. This replacement will give an effective GT strength < GT >'^^ . 
To determine the effective coupling constants, g^, 5(|5i and g^pj^ for a given major shell, one 
takes a sufficiently large number of GT decays that belong to this major shell and whose 
ft values, or equivalently i?(GT)'s, are known experimentally. One then adjusts g^^ 
and g'lpj^ in such a manner that the i?(GT)'s calculated with < GT >'^^ would optimally 
reproduce the chosen set of experimental i?(GT)'s. Once g*^^ , gl^j^ and (7^?^ are determined, 
one can predict all other i?(GT)'s involving the same major shell. 

The EEMO has proven to be highly successful in correlating a great majority of the 
observed magnetic moments. Ml strengths, and Gamow- Teller strengths over a wide range 
in the periodic table |89, 91 1. Furthermore, the effective coupling constants determined in 



EEOM show reasonable agreement with the corresponding quantities obtained theoretically 
for the closed-shell plus/minus one nucleon systems 94 1. 



Two p-shell nuclei, ^^B and ^^C, have recently been considered as useful targets to 
simultaneously monitor the NC and CC reactions for the solar neutrinos. In fact, Raghavan, 
Pakvasa and Brown's proposal |95| to use the ^^B target is already being intensively pursued 
in the BOREX and BOREXINO projects. We describe here how the EEOM (or a closely 
related approach) was used to estimate the cross sections on these targets. 

4.1. 13 c target 

The possible advantage of a ^^C-enriched scintillation counter as a solar- neutrino 



detector was suggested by Arafune et al. |96] and further studied in |^^. These authors 



emphasized that the replacement of ^^C in a scintillator with ^^C not only makes the detector 
much more sensitive to CC reactions, but also allows efficient detection of the NC reactions, 
and that even the natural abundance of ^^C (~1 %) may be large enough to detect NC 
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reactions in a large volume scintillator. The relevant reactions are the CC reaction 

^^C + Ue ^^N + e", (34) 

and the NC reaction 

^^C + u^ "C* + u'. (35) 

For the solar neutrino energy region, only the following final states need to be consid- 
ered For the CC reactions, jJ'^ = 1/2", T = l/2(gnd) >, |3/2-, 1/2(3.51 MeV) >, 
|l/2-,l/2(8.92MeV) >, and 13/2", 1/2(9.48 MeV) > in ^^N, with Q = 2.22 MeV for the 
ground state transition. For the NC reactions, |3/2", 1/2(3.68 MeV) >, |l/2", 1/2(8.86 MeV) >, 
and |3/2-, 1/2(9.90 MeV) > in ^^c. 

The total cross section a for the CC reaction z^g + |^ >— > e~ + \B > is given in 
impulse approximation by 

a = ^^^^PeE,F{ZB,E,) (36) 



vr 

1 



-[| <S|||t|||^> + <B\\\Ta\\\A> f], (37) 



6(2J^ + 1) 

where gA = —1.262. The Fermi matrix element contributes only to a transition between 
isomultiplet members. One often uses the reduced transition strengths -B(F) and i?(GT) 
defined as 

5(F) = J2 \< JBMB\t±\JAMA>\^ = h<B\\\T\\\A>\^, (38) 

Mb, Ma 

B{GT) = J2 \ < JBMB\t±c7\JAMA> \^ = h <B\\\tc7\\\A> \'^. (39) 

Mb, Ma 

The NC cross section is given by 

cr = ^9liElf^\<B\\\ra\\\A>\^ (40) 

For the CC, if the ft value of f3 decay \B >^ |^ > is known, a is given model- 
independently in terms of ft: 
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This applies to ^[^^C N(gnd)]. From log(/t)°^P = 3.667 ± 0.001 for ^3N(gnd) 

^'^C(gnd), we obtain model-independently <a[^^C ^'^N(gnd)]>= 8.12 x 10~^'^cm^. Here 
< a > stands for the cross section averaged over the *B neutrino spectrum. (Due to the 
threshold energy Eth = 2.2 MeV, only the neutrinos from ®B decay are relevant here.) As 
i3]\[(gnd) ^^C(gnd) is a super-allowed transition, one can expect a strong feeding of 
^^N(gnd). For the transitions to the excited states in ^^N, for which no experimental ft 
values are available, one must rely on theoretical estimations. The same applies to the NC 
reactions. 

Pukugita et al. |98| used EEOM to calculate these cross sections. Thus, in eq. (^7|), 



<i?|||T<T|||^> is replaced by 

<B\\\Ta\\\A>^^=<B\\\ffTa + ffATi + ff^T[Y2Xa]^'-^\\\A> . (42) 



Correspondingly, B{GT) will change into B{GTy^ [cf. eq.(|39|)]. {B(F) does not change 
because of charge conservation.) For a model space P, Fukugita et al. used the Cohen- 
Kurath (CK) wave functions [ p9| , which are full-p-shell wave functions. To determine the 
effective couplings, f^, ff^ and /|?4, the following three ground-state transitions were 
used: ^^0{(3^y^N, ^'-^NiP+y^C, and "C(/3+)13B. Since gf^rl is known to be very smah 
@ il> ffA was fixed at a commonly accepted value |9|, H: ff^ = 0.011. The best-fit 



values of the remaining parameters are 

ff = 0.69, ff^ = 0.19. These effective coupling 
constants were used to calculate < S|||Tcr|||A >'^^, i?(GT)°^, and the cross sections for 
the reactions eqs. (|3^), (35). The resulting i?(GT)^''^'s were found to be smaller than 



the corresponding bare shell-model values by a factor ~ 2. According to Fukugita et 
al. [^], in the solar neutrino energy range, one needs to consider only the ground and 
first-excited states in ^^N for the CC reactions, and the first excited state in ^^C for the 
NC reaction The calculated cross sections for the dominant reactions leading to the first 
excited states are: < a[^^C ^3^(3/2" (3.51 MeV))] >= 2.62 x IQ-^^cm^ for CC, and 
<a[^^C 13^(3/2- (3.68 MeV))] >= 1.16 x W'"^^ cm^ for NC. 
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Pukugita et al, comparing ^(GT)^''' for ^3]sj(3/2- [3.51 MeV]) with the corre- 

sponding "empirical" value deduced from the {p, n) reaction, assigned ~30 % uncertainty 
to the calculated a. (The use of {p, n) reactions to determine i?(GT)'s is discussed in chap- 
ter 5.) One could improve Pukugita et aFs calculation significantly by using the results 



of a recent much more elaborate EEOM study |^l|]. Chou et al. |91] have carried out a 
comprehensive analysis of Gamow- Teller P decays for the light (A < 18) nuclei. The model 
space P was taken to be either (Os)^(Op)^-'^ or {Os)^{Op)^'^^'^''' {IsOd)''^'' . The effective 
GT operator was deduced for the Op shell from a least-squares fit to sixteen i3(GT)'^^P's, 
with the results 0.8, ff^ « 0, and /p^i ~ 0. Using the deduced effective operator, 

eighty-three -B(GT)'s were predicted and compared with i?(GT)'^^P's. According to Chou 
et al, their EEOM analysis can reproduce reasonably strong i?(GT)'s belonging to the Op 



shell within ~10 % accuracy. It may be worthwhile to redo the calculation of [98| taking 



into account Chou et a/'s latest result. In this connection one should also consult another 



detailed calculation on the p-shell nuclei due to the Utrecht group [101, 102 1, which utilized 
the full space of Ohuj + 2huj configurations. 

According to Fukugita et al, the ^^C target is expected to have 6.22 x (1 it 0.1) SNU 
from the CC reactions, and 0.68 x (1 it 0.3) SNU from the NC reaction. These correspond 
to 7870 CC events per kiloton yr and 860 NC events per kiloton yr for a ^'^CH2 detector; 
these are both encouragingly large event rates. 

The fact that number of final states that need to be taken into account for the ^'^C 
detector is highly limited implies the possibility of direct calibration using neutrino beams 
from stopped muons. The cross sections averaged over the /x "'"-decay neutrino spectrum is 
estimated to be [Esll 



<a[13C(l/e,e")''N(gnd)]>^+ decay = 2.22 X lO'^^cm^, 
<c7[i3c(i/e,e-)i3j^(^ ; 3.51 MeV)]>^+ decay = 2.27 X 10-41 cm2, 
<f7[i3c(z.e,i^e)'^C(| ; 3.68 MeV)]>^+ decay = 0.56 X IQ-^^cm^ 
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<a[^^C{v^,9^Y^C{- ; 3.68 MeV)]>^+ decay = 0.74 x IQ-^icm^. (43) 

In this estimation, the finite momentum transfer effect is included by multiplying B{F) and 
B{GTY^ with jo{qr). This form factor effect reduces < a >^+ decay by ~ 15%, but it is 
negligible (< 1%) for the solar neutrino energies. These estimates are also relevant to nat- 
ural carbon scintillator experiments, where ^^C contamination may contribute a significant 
background. 



4.2. target — BOREX and BOREXINO experiments ||, [To|, [To|, [To|, |06| 

Borex is a boron-rich liquid scintillation detector that contains approximately 200 t 
of boron in the fiducial volume, and is designed to measure the bolometric solar neutrino 
flux. Borex allows the observation of four basic reactions: 

i/+"B ^ i.' + [iiB*(i?i)^ 7(^x0 +"B], (44) 

^^+iiB ^ e' + [''C*{El)^^{E',)+''C], (45) 

i/ + e" -> u' + e', (46) 

i>e + P — *■ n + (47) 

The NC process involves the three excited states of ^^B at Ex= 2.12, 4.45 and 5.02 MeV, 
while the CC process involves the four states in ^^C at E'^ = 0, 2.0, 4.32 and 4.80 MeV. 
The recent accumulation of solar neutrino experimental data as discussed in chapter 1 has 
enhanced the importance of measuring the solar neutrino flux coming from ^Be, and the 
fact that the ^^B target has an appropriate detection threshold for this purpose further 
increases our interest in the BOREX project. These discussions can be found in []106|| . Here 
we concentrate on the method used by Raghavan, Pakvasa and Brown in their celebrated 
1986 paper [^] to estimate the cross sections for the nuclear processes eqs. (Q), (^5[). 
The CC reaction cross section leading to the ^^C ground state can be directly connected to 



the experimental ft value for the C — > B /3 decay [cf. eq. (41)] . To deduce the cross 



sections for the other final states, eq. ( [3^ ) without the Fermi contribution is used. However, 
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the magnitudes of quenching of the relevant i?(GT)'s are estimated using information on 
the electromagnetic process. That is, by exploiting the facts that the isovector part of 
the Ml operator also involves rcr and that -B(Ml)jgY's between the ground state and the 
relevant excites states in ^^B are known from 7 decay and (e, e') data, the effective B{GT) 
is determined by assuming the proportionality 

B{GT) S(Ml)isv' ^ ^ 

This is in fact the Walecka-Donnelly (W-D) method ||7l|, [7^ discussed in the previous chapter 
rather than the EEOM as defined here. Although the core-polarization effect should quench 
B{GT) and B{M1) equally, the exchange-current effects on i?(GT) can be significantly 



different from those on B{M1). The calculation of the exchange-current effect |93| indicates 
that the use of eq. ( ]48| ) would introduce ~30 % errors. It is therefore in general preferable 
to use EEOM, in which B{GTy^ is deduced using only information on the weak-interaction 
processes. Since the effective GT operator for the p-shell based on the elaborate EEOM is 
now available it seems warranted to derive new estimates of the cross sections relevant 
to BOREX with the use of EEOM. This will probably reduce the uncertainties in the cross 
sections from the 30 % to 10 % level. As far as the classification of calculational methods 
is concerned, we could have discussed ^^B in relation to the W-D method. The reason why 
we discuss it in this chapter is the similarity of considerations that should go into ^'^C and 
"B. 



29 



Chapter 5 Determination of Gamow- Teller strengths through (p, n) reactions 

5.1 Formalism 

The (p, n) reaction has been used extensively to estimate the Gamow- Teller strengths 
-B(GT) that determine the cross sections of neutrino-nucleus reactions of great astrophysi- 
cal importance. The key points here are that the (p, n) reaction, in contrast to (3 decay, is 
not limited by Q-value constraints, and that there exists an approximate "empirical" pro- 
portionality between B{GT) and the forward {p,n) reaction cross section at intermediate 



energies |107, 105, |109| ]. This proportionality was studied in great detail by Tadeucci et al. 
|109| ]. The essential points of their classic paper may be summarized as follows. Consider 
the {p, n) reaction 

p{Ep) + zAn n{En) + z+iAn-i{3), (49) 

where j distinguishes various final states. The transition leading to state j may be of the 
Fermi-, GT-, or mixed type; for the sake of simplicity, however, we first consider the pure 
cases of a-type (a= F or GT), and describe later how to extend the results to the mixed- 
type transition. Let (tJ^''(A, Ep] q, oj) denote the observed differential cross section with uj = 
Ep—En, q = three-momentum transfer, and q = \q\. Since the 0° cross section is particularly 
important, we introduce (T^^^(0°; c<j) = (7^^ [q,Lo)\Q° , where the momentum transfer q for 
0° is related to u as q = ujjvp. For sufficiently high incident energies and for forward 
angles, the direct reaction mechanism should be dominant. Then, for a differential cross 
section o"J'^^'^(A, E'p; q^uS) obtained in a DWIA calculation, we can expect aY^^^isi^^) ~ 
(T°^^(q', Lij). We therefore adopt as a basic ansatz the equality 

afP(0°;u;) = afWiA(o°;^), (50) 

where (tJ*^^'^(0°; w) = crj*^^"^ (g, u;)|o° . In the direct-reaction picture, a set of quantum 
numbers corresponding to the isospin- hipping and/or spin-flipping of the incident nucleon 
is transferred in a single step to the target nucleus. This single-step transfer, which could 
be considered as arising from the exchange of mesons or meson-like effective particles, is 
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completely analogous to what the electroweak current does to a nucleus. We therefore 
expect a^^^^{q,u!) for forward angles to be closely related to the corresponding /? decay 
strength B{a), where a = F or GT depending on j. 

To guide ourselves in quantifying this "close" relationship, we consider a simplified 
DWIA treatment which excludes orbital angular momentum transfer (L = contribution 
only), and in which the initial and final distortion effects are taken into account using the 



eikonal approximation |11C]. This simplified treatment, to be referred to as the eikonal 



DWIA calculation, leads to an explicit proportionality relation |109]: 



af''{q,u;) = K{Ep,co) e^-xA^^^ + p{u;)] | J^j^ exp{-^q^ < rl >)B{a),. (51) 

Although the actual deduction of B{a)j always uses the 0° cross sections, we are considering 
here general forward-angle scatterings, in which q and lv are independent variables. In the 
present context it is useful to treat these two variables separately, since the q- and lv- 
dependence of the cross sections arises from quite different physical origins. The first factor 
in eq. ([5l| ) is kinematic in origin, and its explicit form is 

K{E,,u:) = ^^f. (52) 

The second factor represents the effect of distortion of the initial and final nucleon waves, 
parametrized in terms of the square-well optical potentials Uij = —{Vij+iWij) with radius 
R = r^A^/'^] X = AWir^/vp and p{u}) arises from the difference between the initial and final 
optical potentials. One could calculate p{oo) in the eikonal DWIA, but it will prove more 
useful to determine it empirically. Jq, is, apart from a correction for the antisymmetrization 
of the valence nucleon and target nucleus, the volume integral of the N-N t-matrix in a 
channel corresponding to a: 

Ja = J Va{f')d^r + (Pauli— exchange term), (53) 

where V^(r) is the radial dependence of the relevant t-matrix; t{l,2) = VF(ri2)Ti • T2 or 
t{l, 2) = Vgt(^12) T^i • T2 <?! • (T2- The fourth factor in eq. ( |5l| ) accounts for the momentum- 
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transfer dependence coming from the overall nuclear size effect as well as the g-dependence 
in the L = Bessel transform of Va{r). We factor out the q- and w-dependence in eq. 



as 

afHq,u;) = F^'\q,u;) af^q = 0,u; = 0), (54) 

where 

exp(— -(7^ < >) exp[p(u;) — p{uj = 0)] . (55) 



K{Ep,Uj) 1 2 ^ „2 



K{Ep,uj = 0) 
Similarly, we define 

^DWIA('„ , ,\ 
^'^'''^ = ctDWIA(^ = 0,.; = 0)- 
By definition, F^'^{q = 0, w = 0) = F^^IA^^ ^ ^ g) = 1. The separability of B{a)j in 
^^^^(gja;) [eq. (|5T|)] implies that 

af'^(^, Fp) = (Tf'^(A, Fp; (7 = 0, = 0)/B(a),-, (57) 

is independent of individual transitions j, except for the distinction between the F- and 
GT-types. We therefore write a^^^{A, Ep) instead of aJ''^^{A, Ep). Its explicit form is 

a^^{A,Ep) = KiEp,uj = 0) exp[-xA^/^ + p^uj = 0)] |J„j2. (53) 
We also introduce 

^DWIA(^^ ^ WIA(^^ ^ = 0, ^ = 0)/B{a)j . (59) 

The hatted quantities, a^^^{A, Ep) and (tJ'^^^(^, £'p), which represent {p,n) cross sections 
one would obtain for q = 0, uj = if B^. = 1, are called the unit cross sections. In terms 
of the unit cross section cr^'^{A,Ep), the separability encoded in eq. (51) is expressed as 

^7fk(g,^) = F'''^{q,u:)afHA,Ep)B{a)^. (60) 



32 



Note that, since the way B{a)j enters mto ay^^^{A, Ep) is yet to be determined, the 
definition in eq. ( |59|) is nothing more than a convenient numerical normaUzation; no exphcit 
expression such as eq. ( |58D can be derived therefrom.. 

Guided by eq. (|60|), we try to derive an approximate proportionaUty between B{a)j 
and aY^{A,Ep; q,oj) or, more specifically, between B{a)j and (T'^^^{A,Ep; 0°;w). The first 
step is to check to what extent the eikonal DWIA can reproduce the numerical results of the 
full DWIA calculation. For the purpose of comparing these two approximations concerning 
the reaction mechanism, one can use simplified nuclear wave functions. Tadeucci et al. 
|109| ], employing extreme single-particle wavefunctions, compared the eikonal DWIA and 
the full DWIA over a wide range of A and Ep. Their results indicate that, if the parameters 
appearing in F^''^{q,uj) [eq. (]55|)] are treated as adjustable parameters, F^^^{q,uj) can fit 
F^^^^{q,oj) reasonably well. The effective F factor obtained through this fit is denoted by 
F{q,uj). Thus 

^ exp(-ig2 < ~2 ^^^^^^^ _ ^ q)] _ (g^) 

K [Ep, u; = (J) 6 

where f is the optimal radial parameter, while p{ijj) is the optimal functional form of 
(Quadratic functions suffice to fit the DWIA results.) Furthermore, by fixing the parameters 
in eq. ([5^)] to best-fit (t^'^^^{A, Ep), we can determine the optimized a^^^{A, Ep), which we 



denote by a^XA, Ep). According to Tadeucci et al. |109| , a^fA^Ep) is a smooth function 



of A and Ep and is able to reproduce the average behavior of cF^'^^^{A,Ep) satisfactorily; 
for GT transitions (oj = GT), the scatter around aaj{A,Ep) due to the j-dependence is 
about 7 %. Using these results as well as eqs. ([5^), ([5^), we obtain 



afWiA(,,^) = Ff^'^{q,u)a^'^'HA,Ep)B{a),, 

^ F{q,Lo)aa,{A,Ep)B{a)j. (62) 

Since the various assumptions made in deriving this result are best justified for 0° scattering, 
we specialize ourselves here to the 0° case. Then, invoking the basic assumption eq. (^0|), 
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we arrive at 

afP(0°;u;) ^ F{q,uj)\oo a^^{A, Ep) B{a)j . (63) 

Note that once the type of a transition {aj = F or GT) is specified, the factor F{q, w)|o° cr^j (A, Ep) 
is independent of individual transitions specified by j, apart from the easily calculable kine- 
matical dependence through uj. Thus eq. ( |6^ ) represents an approximate proportionality 
between ctY^{Q°]lo) and B{a)j. The extension of this relation to the case of a mixed tran- 
sition is 

afP(0°;a;) ^ F{q,io)\oo x [aFiA, Ep) B{F)j + <7gt(A ^p) ^(GT),]. (64) 

As a matter of fact, F{q,uj) slightly varies depending on aj= F or GT, through the aj- 
dependence of f in eq. (|6l|). This difference, however, is negligible in practice, and was 
ignored in the above expression. 

In the literature one often encounters another expression which is equivalent to eq. 
(|6^. Dropping the optical potential in what is called here the eikonal DWIA, we consider the 
plane- wave impulse approximation (P WIA) . The unit cross section in PWIA, aJ'^'^^{A, Ep), 
is given simply by ignoring the exponential factor in eq. (|58|): 

aJ'^'^iAEp) = K{Ep,io = 0)\Ja,\^ . (65) 

The ratio 

^ Siiff'J , = eM-^A^' +Pi^ = 0)], (66) 

represents the optical potential effect for the 0° scattering, and is called the distortion factor. 
Eqs. (13), (H) and (|6D lead to 



af^{A,Ep; g = 0,u; = 0) = K{Ep,uj = 0) |J„/S(q)j . (67) 
This is completely equivalent to the 0° value of eq. (|60|). Furthermore, upon optimizing 



the parameters appearing in eq. (l67|) to fit the full DWIA results, eq. (|63) would become 
equivalent to eq. (|63|). 
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The important question is how rehable eq. (p3| ) is. As mentioned above, there is an 
intrinsic uncertainty (~ 7 %) due to the scatter of ay^^^{q,ui) around the smooth function 
aaj{A, Ep). In fact, an even larger ambiguity arises from the optical-model dependence of 
a^^^{q,ui) and aaj{A, Ep). These unit cross sections can vary up to ~ 30 % within rea- 
sonable choices of the optical potential. (The ratio a^^^^{q,u!)/aaj{A, Ep) is not sensitive 
to this change.) Thus, at present, the proportionality coefficient in eq. ( |63| ) cannot be 
calculated with accuracy better than 30 %. 

A useful alternative is to determine the proportionality coefficient empirically using 
transitions for which both aj^^{0°;uj) and B{a)Y^^ are known. For each of these transitions 
one can deduce the empirical unit cross section via 

afP(0°;a;) 



If the proportionality relation eq. ( |63[ ) is a good description of reality, [(T^™'']j's should lie 
close to the smooth curve corresponding to cTq^. (A, £'p). For the GT process, Tadeucci et 
al. [109 1 studied a large number of transitions covering from A = Q to A = 162 for various 



incident energies Ep. Tadeucci et al. refer to the proportionality among transitions with 
a common pair of initial and final nuclides as specific proportionality. The proportional- 
ity within transitions involving different pairs of nuclides is called general proportionality. 
Tadeucci et aVs results indicate that the specific proportionality holds typically to an accu- 
racy of ~ 5%. Thus the use of empirical data helps to reduce ambiguities which at present 
cannot be controlled theoretically. On the other hand, the same authors report that the 
general proportionality may be uncertain by as much as ~ 50 %. In particular, the cases 
involving odd-A targets exhibit the largest deviations from the predicted smooth curve. As 
far as the general proportionality is concerned, therefore, the model consideration described 
here is open to further improvements. The even-odd mass-number effect in the optical 



potential was cited as a possible solution for this problem [109]. 



The above argument suggests as a practical way to use the "proportionality" the 
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following. To infer B{a)j pertaining to a transition in a mass-^ nucleus, find a transition 
j' for which -B(a)^f^ is known from a /3-decay experiment and which belongs either to the 
same mass-^ nucleus or to the closest possible neighboring nucleus. Measure cjJ^^(0°;u;) 
and (T^P(0°;u;). Compute [a^^^lj' from eq. (|6^). Obtain B{a)j using eq. ( |63[ ) together 
with the assumption = This is in fact essentially the method that was used 

by Goodman et al. in their seminal paper |107|, prior to the development of the detailed 



treatment described in [|109 |. 

As a particular case of the specific proportionality Tadeucci et al. also analyzed 
mixed-type transitions [eq. (|6^] leading to the analog states of target nuclei. A very useful 
empirical rule to relate the -B(GT) of an analog transition to the B{F) in the same transition 
was reported in [109]. Note that B{F) is given model-independently by B{F) = N — Z. 



We add that all these results can also be stated using eq. (|67|), which is equivalent 
to eq. (H). 

5.2 Experimental studies 

Systematic experimental studies of the GT-strength in nuclei with the use of intermediate- 



energy {p, n) reactions dramatically started with two seminal articles in 1980 [111, 107]. The 



first paper [111| clearly demonstrated the existence of the giant GT resonance predicted by 



Ikeda, Fujii and Fujita as early as 1963 [112|. In the other paper []107[, Goodman et al. 



|107| ] established the approximate proportionality eq. (33) between i3(GT)'s and the forward 
{p, n) reaction cross section and emphasized that this proportionality could be a powerful 
tool to extract GT strengths for transitions inaccessible to (i decay. Indeed this tool has 



been used extensively in many subsequent works ]108, 109, 113, 114, 115]. Among the main 
objectives of these investigations are: (i) Determining GT-strength distributions over wide 
excitation-energy ranges and comparing them with shell-model predictions as well as with 



the GT-strength sum rule, (see e.g. J116, |115| ]); (ii) Extracting i3(GT)'s for particular 
transitions which play important roles in astrophysical neutrino detections. Although items 
(i) and (ii) are correlated, the latter is of more direct relevance to our discussion, and we 
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mostly concentrate here on (ii). 

A list of neutrino-nucleus reaction cross sections that are relevant to solar neutrino 
detection and that have been determined with the help of (p, n) reaction data includes ^^F, 
3^01 (11^, ^°Ar, ^^Ga ^^Br |T|, |2|, ^^Mo [|l20|, ^^^In, |l2|], ^^^Tl A 



thorough description of specific solar neutrino experiments that use these nuclear targets 
and a critical discussion of uncertainties in the estimates of the relevant cross sections can 
be found in Bahcall and Ulrich and in Bahcall Q. There is no need to repeat what is 
discussed extensively in these references but, in view of great current interest in the results 



of GALLEX 1 13] and SAGE [14, 15 1, it may be useful to remind the reader that the standard 
solar model value of 132l^y SNU for a ''^Ga detector, given in |l| and used in the analyses 
of the GALLEX and SAGE data, already includes rather generous error estimates of the 
B(GT)'s for ^^Ge excited states. [The cross section for the ground-state transition can be 
determined from the /i''^P(^^Ge — > ^^Ga) without using the {p,n) reaction, and hence has 
little uncertainty.] 

We mention here two post-Bahcall-Ulrich developments. One is concerned with the 
z/e + ^^^I ^ e" + ^2^Xe (Q = 0.789 MeV) . (69) 



Originally, Haxton |122(] advocated the usefulness of an iodine detector for testing whether 
the solar neutrino flux is constant as a function of time. Our interest in the iodine detector 
has recently received an additional boost from the fact that it is sensitive to the ''Be neu- 
trinos as well as to the ^Be neutrinos. The results of the GALLEX and SAGE experiments 



p3| , 14, |15[ have further enhanced the importance of examining the energy dependence of 
solar neutrino deficits with the use of nuclear targets with different reaction Q- values. The 
^^''l target is a promising candidate to cover the region between the ''^Ga experiment and 
the water-Cerenkov-counter and ^^Cl experiments. The main issue with ^^'^I experiments is 



how to estimate the cross section of eq. (p9D reliably, and this problem is currently studied 



with great intensity. ]123]. For an update of (p, n)-reaction determination of the i3(GT)'s 
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for ^'^'^1 as well as for^^Ga, ^^Br, see e.g. Sugarbaker [124|. For the current status of cali- 



brating the ^^^I(i^e, e~) ^2^Xe cross section with the LAMPF muon-decay neutrinos, see e.g. 



Wildenhain [125]; for the nuclear structure consideration of the relevant i?(GT)'s, see Engel 



et al. Ill] and Haxton fTl^ 



Another topic that has recently been attracting much attention will be discussed in 
the next section. 

5.3 The low-energy Vf, + ^^Cl — > e~ Ar reaction 
As is well known, the reaction 

Ve + ^^Cl ^ e" Ar, (70) 

was the first to be used for detecting the solar neutrinos. It is also well known that the 
determination of the relevant cross sections presents a unique case. That is, there exists 
decay 

37Cai7^?^Ki8 + e+ + i/e, (71) 
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which is the mirror process of 



?7Cl20^?lAri9 + e- + i>e, (72) 



and all of the '^''Ar states that are relevant to the solar neutrino reaction eq. (|7^) have their 
counterparts in the final ^'^K states of the ^''Ca /3^-decay eq. (^), since this decay has a 
very large energy release ((5ec=11-64 MeV). Therefore, using the experimental /f-values for 
individual final states in eq. (|7^), one can determine model- independently low-energy cross 



sections for eq. ( |7(]| ) [c/. eq. (|4l|)]. The mirror-asymmetry problem J128] can in principle 
affect this procedure to some extent. The mirror asymmetry parameter 5 is defined by 
6 = {ft)^/{ft^) — 1 for a pair of mirror and (3~ decays. With no charge-asymmetric 
effects in nuclear structure, and in the absence of second-class currents (SCO) in the weak- 
interaction current, 5 should be identically zero. Experimentally, 5^^^ ^ 0.1 for even-^ 



systems, and (5™p can be as large as ~0.2 for odd-^ nuclei [128| . (It is commonly accepted 
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that the finite values of J'^^p are caused by the "trivial" nuclear charge-asymmetry effects 
rather than by the SCC; no explicit calculations that can explain the large (5°^p in the odd- A 
cases have ever been presented, however [128, |129| .) The 6 effect should in principle exist 



in the A = 37 case as well. Fortunately, however, the solar-neutrino reaction on '^''Cl is 
dominated by the feeding of the analog-state {Ex= 4.98 MeV) whose strength is primarily 
governed by the well-known B{F) = N — Z, a situation that renders the estimation of 
a['^'^Cl{ue, e~)^'^ Ar] significantly less sensitive to the mirror-asymmetry problem than the 
general magnitude of 5^^^ may suggest. As one goes beyond the solar neutrino energies, 
however, the contributions of the higher-lying states become progressively more important 
and hence the 6 problem may acquire a practical significance. 

Another problem whose importance increases as becomes larger than the solar 
neutrino energies is the discrepancy between the -B(GT)'s deduced from the ^"^Ca /3 decay 
| 13C ] and those obtained from the {p, n) reaction 117| . In discussing this problem, it is useful 



to consult the shell-model prediction |p9| obtained with the use of the empirical transition 



operator method (cf. chapter 4) based on the Wildenthal wavefunctions |131| . Now, while 
^"^Ca ^"^K + Z?"^ data |130| give i?(GT)'s in fair agreement with the shell-model values 



for < 5.7 MeV, they yield too smah S(GT)'s for the levels in the E^^ = 6.5 - 8.5 MeV 



range. According to Adelberger and Haxton |132], this discrepancy arises because Sextro 
et al. ]130| ] assumed that the delayed proton emission from '^^K leads only to the 0^ ground 
state of ^''Ar, thus ignoring the expected substantial feeding of the 2"'' (3.827 MeV) state; if 
the decay to this 2"*" state is properly taken into account, the resulting i?(GT)'s are rather 



close to the shell-model prediction. The {p, n) data |117] indicate a strong (summed) B{GT) 
in the E^ = 6.5 — 8.5 MeV region in agreement with the shell-model prediction. The overall 
agreement between the -B(GT)'s obtained from the {p,n) reaction [117| and the shell-model 



values [^] is reasonable except for levels in the vicinity of the analog state [E^ = 4.98 
MeV), for which the '^''Cl(p, n)^''Ar reaction data give too small S(GT)'s. Adelberger and 



Haxton [132| suggested that this depletion of the GT strengths was caused by an incorrect 
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subtraction procedure of the strong Fermi contribution. 



Garcia et al. [133] recently carried out a high-resolution, low-background experiment 
on '^''Ca decay, drastically improving the time-honored work of Sextro, Gough and Cerny 



|13C]. For Ex < 5.5 MeV, the -B(GT)'s obtained by Garcia are in fair agreement with those 
of [13C], except that Garcia et al. see no evidence for feeding of a 4.679 Mev level in ^''K. 



On the other hand, Garcia et al. find much more strength at Ex > 5.5 MeV, and the 
integrated B{GT) is roughly twice that of Sextro et al. Furthermore, Garcia et al. identify 



~50 % more GT strength below Ex=S.Q MeV than was inferred from the {p,n) data [117|. 
According to Adelberger et al. [ [134| , the integrated B{GTY^'^ up to Ex=^ MeV does not 
require the renormalization of qa (see also [|115[1 ), and the quenching of qa indicated by 



shell- model analyses |90, 116| of the (3 decay data may in fact suggest that the shell-model 



calculations tend to shift the i3(GT)'s to the lower excited levels. Since the /3-decay 

data concerns only low-lying states, the shell-model fit of the observed /3-decay strengths 
leads to the apparent quenching of qa- 



Critical examinations of various aspects of Garcia et aVs analyses |132, 134 1 have 



been made in |135| , 136, 137|. The comparison of the latest ^"^Ca decay result with the 



(p, n) data points out among other things the following two problems: (i) On the (p, n) 
side, the determination of the B{GT) in the isobaric-analog state region {Ex ~ 5 MeV) 
is quite delicate because of the dominant Fermi transition and the limited experimental 



energy resolution [|135|| ; (ii) On the /3-decay side, the assumption used in 133 | that all the 
^"^K daughter states above the proton threshold should decay by proton emission is not 
tenable |3|, |l3|]. In fact, the 3.24 MeV level in has F^/Fp « 40 [|4|], and this can 



change significantly Garcia et aVs analysis of the delayed-proton emission. 

A detailed comparison between the updated values of B{GTy^^ and the improved 



shell-model prediction was made by Brown |139| . The results of a recent re- measurement 
1 135] of the ^''Cl(j», n)'^'^Ar are expected to further clarify the issue. 

It is to be remarked that, although from a nuclear-structure point of view it is im- 
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portant and interesting to settle the existing quantitative difference between tlie ^^Ca j3~^ 
decay result and the {p, n) data, the influence of this discrepancy on the solar neutrino cross 



section is rather modest; for instance, the -B(GT)'s deduced by Garcia et al. |133| implies 
that (T[^'^Cl(i^e, e~)^'^Ar] for the neutrinos would be enhanced by 6% over the literature 
value. 



41 



Chapter 6 Inclusive neutrino-nucleus reactions - Fermi gas model 

6.1 Overview 

Our interest in intermediate- and high-energy neutrino-nucleus reactions was greatly 
enhanced by an intriguing issue raised by Davis 0. According to Davis, the '^''Ar production 
rate in the CI experiment has a significant time structure which may be correlated with the 
occurrence of large solar fiares A possible explanation is that, when protons accelerated 
by the flare activity collide with the solar atmosphere, they may produce a sufficient fluence 
of medium- and high-energy neutrinos (to be called the solar-flare neutrinos) to cause signif- 
icant excess counts of '^''Ar. The typical energies of the solar-flare neutrinos are Ey a. few 
hundred MeV, and their high energy tail can even extend to the GeV region. The original 
suggestion ||^] was based on a crude estimate of the z/-^''Cl cross section at Ey =100 MeV, 
but quantitative considerations would require more elaborate estimates. Davis's argument 
based on the Homestake experiment was followed by the suggestion @, |4l|, |l4|, |l5l that 



large water Cerenkov detectors such as the Kamiokande and 1MB facilities can furnish a 
powerful tool to measure the fluence of solar-flare neutrinos. The point is that because 
of the high energies of the flare-neutrinos the ^^O nuclei that are abundantly contained 
in a water detector become extremely efficient targets for detecting neutrinos via inclusive 



reactions [25|. A similar situation exists for the ^^C nuclei in scintillation counters [ 141 |. 
These arguments call for reasonably reliable estimates of neutrino-nucleus reaction cross 
sections for Ei^ typically up to a few hundred MeV. Two different problems are involved 
here. For radio-chemical experiments that count only those events that lead to particle- 
bound states, we must calculate semi-inclusive cross sections, whereas water Cerenkov or 
liquid scintillation experiments involve the usual inclusive reaction cross sections. Relegat- 
ing the discussion of the semi-inclusive reactions to chapter 7, we survey here recent works 
concerning the inclusive reactions. 

The inclusive neutrino-nucleus reactions are also important in relation to the atmo- 
spheric neutrinos. When primary cosmic rays strike the atmosphere, they produce pions and 
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kaons that subsequently decay into muons and muon-neutrinos, and much less frequently 
electrons and electron-neutrinos. The muons further decay into electron-neutrinos and 
muon-neutrinos. The neutrinos of this origin, called the atmospheric neutrinos, have been 
observed in the range of a few hundred MeV to several GeV with the large water Cerenkov 



detectors [143, 144]. The interactions of the atmospheric neutrinos in water are important 
for two main reasons. One is that these interactions constitute the principal background for 
nucleon decay experiments and extraterrestrial neutrino searches. The second reason is an 
important current issue concerning the atmospheric v^/ve ratio. The above-described ori- 
gin of the atmospheric neutrinos implies that there should be roughly two muon-neutrinos 
for each electron-neutrino. (Here "neutrinos" mean both neutrinos and antineutrinos.) 
The experimental data |143| , |144| ], however, indicates that the muon-neutrinos are signif- 
icantly less abundant than expected. The cross sections for ^^0(2^^, £~)X(anything) and 
^^0(Pf , (anything) feature importantly in analyzing these data. 

6.2 Calculational methods 

The study of inclusive nuclear responses to medium- or high-energy electroweak 
probes has a vast accumulation of literature |145| , [n| , [7^ , 14C, 147, 148 , 149 , 151 1. The 



closure approximation | 146 , 147 1 and the quasi- free scattering approach [jTll , 148 , 149 , 151 1 



are the most commonly used methods. Kuramoto et al. [154, 155] have recently made a de- 
tailed comparison of these two formalisms and recommend that the closure approximation 
not be used for quantitative estimations. We therefore concentrate here on the quasi-free 
scattering approach. 

We are concerned with the inclusive reactions 

1^1+ A + X (anything), (73) 

i^e +A ^+ + X (anything). (74) 



In the quasi-free scattering (QFS) approach 145 , 71] the transition amplitudes for the 



above processes are obtained by summing the transition amplitudes for neutrino-nucleon 
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scattering over A nucleons characterized by the initial and final nuclear wavef unctions. 



In the so-called relativistic treatment [148] one uses the original single nucleon amplitude 
given in terms of the Dirac spinors In the non-relativistic approximation, one uses 

the non-relativistic (NR) reduction (the Foldy-Wouthuysen transformation) of the single- 



particle transition amplitude 



In the NR approximation the cross sections for eqs. ([73|), (|7J) can be factorized into 
the single-nucleon cross section and the nuclear response function which accounts for the 



nucleon distribution inside nucleus and the Pauli blocking effect. Bugaev et al. |149| used 
the NR form truncated at first order in \q\/M, where q is the three-momentum transfer, 
and calculated the inclusive cross sections on ^^C, ^^O, ^''Cl, ^^Ga and ^^Br, for ^ 300 
MeV. The nuclear wavefunctions were taken to be shell model wavefunctions in suitably 
chosen Woods-Saxon potential. Both bound and unbound orbits were allowed as the final- 



state orbits. Gaisser and O'Connell |151] used the Fermi gas wavefunctions and calculated 
the inclusive cross sections for Ei, = 50 MeV ~ 2 GeV. To the extent that the Fermi 
gas model parameters {pf{p), Pri'iT'), and the average nucleon separation energy Eb) can 
characterize actual nuclei, the results can be applied to any nuclear targets. Above the 
pion production threshold, X in eqs. ([73|), (^) should include meson production processes, 



but they were not considered in |151]. (For neutrino- nucleus reactions leading to meson 



productions, see |152, p,53[ .) According to Gaisser and O'Connell |151| , (i) the relativistic 



and NR treatments give virtually identical results, (ii) in the energy region where their and 
Bugaev et aVs calculations overlap, the results are similar, and (ii) the results are accurate 
to energies as low as 50 MeV for electrons, and to threshold for muons. Thus, although 



Gaisser and O'Connell's calculation |151] has been playing an extremely important role in 



analyzing the f^/ve ratio problem |143, 144 1, the measure of possible errors in the results 



is not readily available from | 151 ] 



Kuramoto, Fukugita, Kohyama and Kubodera [155] made a detailed study of the 



i/(z^)-^^0 inclusive reactions for Ei, ^ 300 MeV, using the relativistic Fermi gas (RFG) 
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model of Smith and Moniz [148| as well as an improved non-relativistic Fermi gas (NRFG) 



model that incorporates the Foldy-Wouthuysen transformation up to third order in \q\/M. 
One of the motivations for Kuramoto et al. to go beyond the lowest-order approximation of 



1 151] is that in some kinematical regions the differential cross sections given in |151] become 
negative, a rather disturbing aspect. Kuramoto et aVs NRFG treatment is found to be free 
from this problem and in reasonable agreement with the RFC model results. The NRFG 
cross sections systematically tend to be somewhat smaller than the RFG results and the 
difference, which grows with E^, reaches ^10 % at = 300 MeV. The agreement between 



the results of Kuramoto et al. and those of Bugaev et al. [149] is also at the 10 % level. 

For lower incident energies, there were attempts to add up explicitly the contributions 
of individual final nuclear states. Langworthy, Lamers and Uberall (LLU) |156| ] considered 
eighteen final states for z/^+^^O £~ (anything), and used the generalized Helm model 
to parametrize individual transition strengths, which were determined semi-empirically with 
the use of the electron scattering data. Haxton ]^] carried out a summation of all possible 
final ^^F states contained in a given truncated shell-model space; all configurations up to 
2huj excitations were included for positive-parity final states and those up to Ihuj excitations 
for negative-parity final states. As in the empirical effective operator method (cf. chapter 
4) the shell-model transition rates were rescaled by adjustable multiplicative factors (up 
to ~ 70% correction) so that the empirical -B(GT) values and the electron scattering data 
were reproduced. As discussed in Kuramoto et al. ]]155(] , fTHaxton agrees very well with ctrfg 
for El, ^ 60 MeV but starts to deviate appreciably beyond E^ 60 MeV; (Tllu exhibits a 
similar tendency, but its deviation from (Trfg beyond Ei, ~ 60 MeV is more pronounced. 
Attributing these features to the fact that the LLU |156| ] took into account only eighteen 
final nuclear states whereas Haxton included much more, Kuramoto et al. argue that the 
difference between fiHaxton and (Trfg for Eiy ^ 60 MeV is due to the limited size of the 
shell-model space used by Haxton. A later work of Haxton and Johnson ]|157[] drastically 
extended the shell-model space and calculated the /3-decay and /x-capture. It would be 
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interesting to see what cross section this extended calculation will give for the inclusive 
reactions, eqs. (|73|), (ff^). 

Although it is gratifying that the RGF results smoothly join those of the explicit enu- 
meration method around ~ 60 MeV, all these methods share the basic assumption that 
the quasi-free scattering picture is valid. The recent data on the transverse and longitudinal 
nuclear responses, Rt and Rl, in quasi-elastic (e, e') reactions indicate significant deviations 
from the Fermi-gas prediction even in the quasi-elastic peak region; the discrepancy can be 
as large as ~ 20% for ^^C and ~ 50% for ^'^Ca. Thus, the original remarkable success 



1 155, 15£, 16C] of the Fermi-gas model in predicting the total (transverse -|- longitudinal) 
response needs to be taken with caution. There have been many attempts to explain Rt 
and Rl, but a consistent picture is yet to be found [161]. The fact that we are dealing here 



with the total nuclear response is a favorable factor, but the fact that we are calculating 
the angle-integrated cross sections that include the non-quasi-elastic region may spoil to 
some extent the quasi-free scattering picture. Taking the mass-number dependence of the 
above-mentioned deviation from the Fermi gas model, Kuramoto et al. assigned ~ 30% 
errors to their estimation of the total cross sections for U£ O ^ £^ + X(anything) and 
ue 0^£+ + X (anything). 

Singh and Oset |162(| have recently made an extremely detailed study for ^^^C and^^O, 
using another version of the Fermi gas model, the local Fermi gas model, which was tested in 
several other processes. Furthermore, these authors considered the nuclear-medium quench- 
ing of the weak-interaction effective coupling constants and the final-state Coulomb distor- 
tion effect. Comparison of Singh and Oset's results with those of Kuramoto et al. also 
indicates that the Fermi gas model is probably reliable up to the ~ 30 % level. 

In actual experiments using a water Cerenkov detector, i/£ and Pi can also scatter 
from electrons. It is however easy to separate these events from v (p) scattering off nuclei, 
since the former has a very sharp forward peak within 9 « (me/E^y^'^. Also, if the incident 
flux consists of a roughly equal mixture of Ve and v'e, the event rates of the reactions off 
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oxygen nuclei will compete with that of z^e + p ^ e"*" + n, and it is important to separate 



one from the other. This can also be done using the angular distribution | 155 ]. 

The significance of the improved estimates of the cross sections for ^^0(f£, £~)X(anythin 
and ^^0{i>£,£~^)X {any thing) for the solar-flare neutrino problem will be discussed at the end 
of the next chapter. As far as the atmospheric neutrino problem is concerned, the question 
is to what extent the ambiguity in the Fermi gas model used in interpreting the data would 
affect the conclusion. Although one generally expects that a substantial part of the ambi- 



guity will disappear upon taking the I'^/i'e ratio |16S], a more quantitative consideration 
seems warranted. 
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Chapter 7 Semi-inclusive neutrino-nucleus reactions 

7.1. 

When high-energy neutrinos hit a nuclear target used in a radiocliemical experiment, 
many nuclear levels can of course be excited but the experiment by design registers only 
the reactions leading to particle-stable states. If the number of final bound states is small 
enough we can hope to calculate the contributions of the individual levels and sum them 
up. (This might be the case for very light nuclei.) As Ei, gets higher, however, the possible 
final bound states can be too numerous for the explicit enumeration of the individual levels 
to be practical. We call this type of process the semi-inclusive reaction. Semi-inclusive 
reactions cannot be treated in the ordinary inclusive-reaction formalisms since particle- 
unbound states must be excluded from the final states. A practically interesting example 
of semi-inclusive reactions was encountered when Davis 0] considered the reaction between 
the solar-flare neutrinos and the ^'^Cl target: 

+ 37ci ^ r + 3^A = e or (75) 

where the final ^''A stands for all particle-stable states. Davis's original argument used a 
rather crude estimate of this semi-inclusive cross section, which motivated Kuramoto et al 



|154| , 155 1 to make a more quantitative estimation. Although Kuramoto et aVs treatment is 
highly specific to the A=37 system, we briefly describe that work here as a possibly useful 
example, since semi-inclusive reactions are likely to be relevant to other radiochemical 
experiments as well. 

The general expression for cross section a for the reaction z/^ + |i >— > £~ + |/ > is 

a = (27r)4 Yl I dhiKEi + Ef-E,-Ei)\< ({p,); f\H,s\u{p,y, i > \\ (76) 
/ 

The final nuclear states \f > should lie below the neutron threshold energy, E^ = 8.8 
MeV. In impulse approximation to be used here, ifeflf is the sum of the single-particle 
weak-interaction operator. The precision of the calculation described here does not warrant 
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taking account of the exchange current effects and other fine details. In this approximation 

a = —cos^ecy2peEeF{Z,Ee)l f' d(cos0)M^, (77) 

where Eg = — AEji {AEfi = Ef — Ei) is the outgoing lepton energy. The squared 
nuclear transition matrix elements Mp is given by: 

Mp = Mf\ <f\l\i>\'' + MGo\ \ < > P + Mg2A (78) 



where 



1 ^ 



<fW>? = r E \ <f-,Jf,Mf\Y,U{k)e'1-''^\i-J,,Mi>\'' (79) 

^-^^ + ^ Mf,AU k=l 



k 



<mi>? ^ ^j— Y E \<f;Jf.Mf\Y,t+{k)a{k)e'1-'->=\i-J,,M,>\''m 



^■^^ + - Mf,M, k=l 

= VrT^T.\< Jf\\j:^+ik)Mqn)[UH) X ^(A;)]W||J, > |^ (82) 

^■^^ ^ £,K k 

and A represents the interference between the < > and < f\&\i > terms. The coef- 
ficients, Mp, Mgo and Mg2, contain the nucleon form factors as well as the kinematical 



quantities [155|. Now, of the transition operators appearing in eqs. (|80|), (|82D, the operators 
that satisfy the selection rule: "AJ = or 1, Avr = no" are classified as allowed-type op- 
erators. All the other operators are grouped together as forbidden-type. (This classification 
slightly differs from what is customarily used in /? decay theory.) In the present case with 
\i> = |^''Cl(gnd); |^ >, the spin-parity of particle-bound states \ f> that can be reached via 
allowed-type operators is or The assembly of these states will be referred to as 

group A. Meanwhile, particle-bound states that can be fed only through the for bidden- type 
operators constitute group B, which contains negative-parity states as well as higher-spin 
positive-parity states. 

For those levels of group A that predominantly belong to sd-shell configurations, 
reasonably realistic wave functions are available from the work of Wildenthal |131|. These 
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wave functions are obtained by diagonalizing full sd-shell configurations and, with the use 
of the empirical effective operator method, the Wildenthal wave functions generally lead to 



a satisfactory agreement with the observed i3(GT)'s |131, The result of this detailed 
rescaling of the strengths can be simulated reasonably well by replacing the free-nucleon 
value gA = —1.262 with an effective value g'^ « —(0.7 ~ 1.1) [132|. At the time of Kuramoto 



et a/.'s work, both the B{GT) values determined from the ^'^Cl{p, n)^'^Ai reaction and those 
deduced from the '^'^Ca ^"^K + (3~^ decay had some irregularities either among themselves 



or as compared with the shell-model predictions. Kuramoto et al. |155| used as the reference 
GT strengths the average of the -B(GT) values determined from these two methods, except 
that, for the Ex > 6.5 MeV region, they adopted the (p,n) data in view of the Adelberger and 
Haxton's argument |132| and, in the ~ 5 MeV region, they used the i?(GT)'s obtained 
from the ^^Ca /3-decay with their strengths slightly reduced according to [|132||. Comparison 



of this reference B{GT) with the shell model prediction indicates that the shell model value 
(with gf « -1.00) slightly overestimates B{GT) for E^ < 4MeV (by a factor of 1.38 on 
the average) and slightly underestimate it for E^ > 4 MeV (by a factor of 0.78). In order 
to reproduce B{GT) for E^ < 4 MeV, one needs to use g*^ = —0.85. Kuramoto et al. 
assumed that this coupling constant could be used universally for the whole range of Ex- 
This implies that the average shell- model value of B(GT) for E^ > 4 MeV is smaller than 
the experimental value by a factor of 0.78 x 0.72 = 0.56. This deficit was interpreted as 
representing the contribution of configurations outside the sd-shell. Using this information 
on B{GT), Kuramoto et al. calculated a a as 

£;^<4MeV 4<£;^<8.8MeV 

where all the axial-vector contributions were calculated for the Wildenthal wave functions 
using g*^ = —0.85, and a = = 1.77 was introduced to "recover" missing contributions 
of the outside-sd-shell contributions. 

The levels belonging to group B may be divided into two subgroups: group Bgd con- 
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sisting of positive-parity states with J/ > | which belong to the sd-shell, and group B' that 
comprises all others. The contribution of group B^^, denoted by a'^, is estimated using 
the Wildenthal wave functions and g*^ = —0.85. The contribution of group B', denoted 
by aB', is assessed by using as' = (t^°'^^^'^Nb' , where a^°'^^^'^ is an estimate of the typical 
cross section to excite a level of group B', and N^' the estimated number of levels belong- 
ing to group B'. The 0-*^°'^^''^ was estimated by taking the average of a^^ and ((T2f)forbid) 
where (u^'^f) forbid is defined through the decomposition a'^/, = ((T2^)aiiowed + (<7i'') forbid, with 
('7A')aiiowcd [(<72^)forbid] being the contribution of the allowed-type [forbidden-type] opera- 
tors. The assessment of N^' is a delicate job, but referring to the available information 
1 164] on the level structure in ^''Ar, Kuramoto et al. estimated Nb' ~ 28 ~ 34. 
The total contribution is obtained as 

a = aA + CTB = CTA + crl^ + ctb'- (84) 



Using this expression, Kuramoto et al. ]155] calculated the cross section for the semi 



inclusive process eq. ([75|). The cross section rises up to 150 ~ 200 MeV and then levels off. 
This is a generally expected behavior for a semi-inclusive reaction feeding bound states, since 
the growth of the phase space factor is counterbalanced by a suppression due to the nuclear- 
size effect. Kuramoto et al. estimated the error of their results to be 20 % at 50 MeV, 40-50 
% at 100 MeV, and 50-60 % or more at 150 MeV. The cross section for > 150 MeV may 
be subject to a factor 2 uncertainty. For the use in calibration experiments with a muon 
beam, the cross section averaged over the neutrino-spectrum from the stopped muon decay 



is relevant. Kuramoto et al. yi55|] give o-('^''Cl)^+_dccay = (9.4 ± 1.4) x 10 "^-^cm^. This 



value should be compared with 9.90 x 10 cm given in 1165], and with 7.6 x 10 cm 



Although the method used by Kuramoto et al. is far from systematic, it is at present 
difficult to do better. It is hoped that the progress of main-stream shell model calculations 
will eventually reach the stage wherein the adopted model space is sufficiently large to 

51 



accommodate all important final states of semi-inclusive neutrino-nucleus reactions, at least 
up to sd-shell nuclei. These shell-model model wavefunctions plus a judicious use of the 
empirical effective operator method (chapter 5) will probably allow us to treat the semi- 
inclusive reactions in a more systematic way. For nuclear targets beyond the sd-shell, the 
possibility of such a systematic treatment seems rather remote. 

7.2. Consequences for the solar-flare problem 



Kuramoto et al. |155] used their estimate of the semi- inclusive reaction eq. ( |75[) 
to deduce the solar-flare neutrino flux corresponding to the fifty excess ^'^Ar events in 
the 615 ton C2CI4 Homestake detector; this excess was assumed by Davis [Q] to be of 
the solar-flare origin. For the flare-neutrino spectrum, a soft neutrino spectrum of the 
form (j){E^) = (po exp{—E^/EQ), with Eq = 50 MeV, as well as a monoenergetic form 
with various plausible energies was used. Based on the solar-flare neutrino flux deduced 
therefrom and their estimate (chapter 5) of the cross section for ^^©(t'f, £^)X(anything), 
Kuramoto et al. calculated the number of neutrino events to be expected in the Kamiokande 
detector during the solar flare periods and concluded that the observed counting rate of the 



neutrino events [143] fell far below the prediction. This corroborates the earlier conclusion 
of the Kamiokande group [|142(| , that the time-integrated "solar-flare fg" flux is less than 
3.7 X 10^cm~^ per flare for E^ = 100 MeV, which is 2000 times smaller than the value 
required to account for the excess "signal" observed in the ^"^Cl experiment. Thus, the 
role of the solar-flare neutrinos in giving rise to the time structure in the Homestake solar- 
neutrino data must be much smaller than described in IItII . For further detailed studies on the 



statistical significance of the time structure in the solar- neutrino data, see |16£, 167, 168 1. 
The untenability of Davis's original suggestion does not diminish the significance of using 
water Cerenkov counters, like the Kamiokande or 1MB, to monitor intermediate energy 
neutrinos associated with the solar activity. This will allow us to explore the active Sun 
with a new, unconventional means, and also to have a better understanding of the particle 
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acceleration mechanism around the solar sphere. It is to be emphasized that the sensitivity 
of the Kamiokande and 1MB water detectors reaches the flux level of an optimistic estimate 



for a large flare or for more constant activity during the solar maximum period |16£, |17C 
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Chapter 8 Additional remarks 

In the foregoing chapters we have discussed various nuclear physics aspects that 
go into the estimation of neutrino- nucleus reactions of astrophysical interest. Neutrino- 
nucleus reactions can also be useful in studying the structure and symmetry properties of 
the fundamental coupling of the weak current with matter. In some of these experiments 
nuclei serve just as an assembly of A nucleons, whereas in many others one takes advantage 
of a wide choice of selection rules pertaining to specific nuclear targets. 

An example of the former is the recent determination of a stringent upper limit to the 



second-class current by Ahrens et al. |171]. These authors measured quasielastic reaction 
Ufj^p fi^n occurring in the scintillation target detectors for up to 1.0 (GeV/c)^. For 
the second-class axial-vector current strength defined by 

<n{pf)\Al^^\p{p,)>= ■i^u{pf)^5'J,.uq^u{pi) (85) 

with q = Pf — Pi, and for the dipole parametrization of the g^-dependence frio^) = 
r//A(0)/(l +q^/M^f, Ahrens et al. deduced \r]\ < 0.8, 0.23, 0.12 for Mr(GeV/c2) = 
0.5, 0.9, 1.5, respectively. These upper limits are substantially lower than the value ob- 



tained from analyzing the totality of the existing (3 decay data. [172, 66]. To the extent 
that the incoherent sum of the individual quasielastic scattering off a nucleon provides an 
accurate description of the process, no nuclear dynamics is involved here. Although this 
picture is likely to be valid, it is probably worthwhile to check its validity in a quantitative 
manner in view of the breakdown (up to 30%) of the quasi-free scattering assumption in 
the longitudinal and transverse response functions in (e, e') processes. 

There have been numerous attempts at exploiting various selection rules available 
from specific nuclear quantum numbers to study unconventional pieces of the weak current. 
Although it is beyond the scope of this review to discuss them in detail, we mention here as 
an example a possible consequence of the strangeness content of the nucleon for neutrino- 
nucleus reactions. Although this effect will not change drastically the interpretation of 
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astrophysical neutrino-nucleus reactions, it has its own fundamental importance. 

The strange quark matrix elements of the nucleon are one of the central issues in 
contemporary particle-nuclear physics. The measurement of polarized deep-inelastic muon 
scattering on the proton [|173|| , combined with an analysis of hyperon /3-decay [174], suggests 



that there is a large contribution of the strange sea quarks (ss) to the proton spin. The 
pion- nucleon sigma term 'SjrN also indicates a significant amplitude of < A^|ss|A^>, although 



the latest estimate [17£] of T^t^n favors somewhat smaller ss content than in |175]. Beise 
and McKeown |177] ] have discussed parity violation measurements of polarized electron- 
proton scattering to study strange quark matrix elements of the vector current. If such a 
large ss component exists in the nucleon, one should expect a sizable contribution also in 
the neutrino-nucleus reactions. In fact, neutrino neutral-current reactions offer the best 
means to study the strange quark matrix element of the axial current, which complement 
the electron-scattering experiments that are sensitive only to the vector current. 

In the presence of the strange quarks, the hadronic neutral current appearing in 
eq. (p^) would become J^^, where 

= + ' (86) 

with the isovector and isoscalar pieces defined by 

Jj=^ = il-2sm^ew)V;^ + Al (87) 
Jj=° = -2 sin2 Ow + + . (88) 

The standard forms for the effective isovector currents in the nucleon sector are 

= gAUNnfii5^UN , (90) 

where Fi and F2 are isovector form factors, and 5a(0) = —1.26. Similarly, the isoscalar 



and the strange quark currents |177| are expressed as 
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-^if' = ^GfflA,.7^7,4"A', (93) 

where and -F2 isoscalar form factors, and Ff, and Gf denote the strange quark 
form factors of the nucleon for the vector current and for the axial vector current, respec- 
tively. Attempts to extract the strange quark axial form factor Gf (0) have been made using 



both EMC [173, 174] and u-p scattering data |178, 179]. The extracted form factors are 



-0.38 ±0.11 from EMC 

cm = { (94) 
-0.30 lb 0.16 from v-p. 

However, it should be emphasized that the quoted numbers are not yet well established. 



because of several assumptions made in the analyses ]17S], so that the quoted errors are the 
experimental ones only. An interesting question is to what extent neutrino-nucleus reactions 
can sharpen our conclusion on the strange-quark neutral current form factors. Henley et al. 
|180| ] have investigated the asymmetry between the elastic cross-sections for u + d ^ u + d 
and P + d — > P + d. In impulse approximation, the difference in cross sections is given by 



i2 



= -^^^v{1-C)GI{f? + F^)f^, (96) 

where rj = q^/Am\i^, C = /^"^NEy, with q the momentum transfer. Ey is the neutrino 
(anti-neutrino) energy, and F£) is the deuteron form factor. According to Henley et a/.'s 
calculation the asymmetry R defined by 

^ {da/dq^){vd vd) - {da/dq^){i>d ud) 

~ {da/dq^){ud ^ i^d) + {da/dq^){ud ud) ^ ' 



could be as large as 15 % for large q^ (> 0.5 GeV^/c^). Suzuki, Kohyama and Yazaki]181] 
made a detailed study of low-energy u-^'^C and u-^^B reactions to extract information on 



Gf . Bernabeu et al. [182] studied the effect of the strange quark matrix element for the 
neutral-current reaction 

+ 1/ + i^B* . (98) 
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This reaction is particularly interesting, because it is directly related to BOREX or BOREX- 
INO experiment (9^. It is to be noted that, unlike the ud reactions, the isovector and 
isoscalar currents contribute coherently to the cross section here. A similar situation should 
exist for the i/^^C neutral-current reaction discussed in chapter 4. 

The above considerations were all based on the impulse approximation. Very recently. 



Hjorth- Jensen, Kirchbach, Riska and Tsushima [183| have pointed out that the exchange 
currents can make the effective isoscalar axial current density of a nucleon in a nucleus 
considerably different from the corresponding free nucleon value. The isoscalar axial charge 
coupling (time component) can be enhanced by 40-50 % due to the short-range component 
of the nucleon-nucleon interaction, whereas the isoscalar axial current coupling (space com- 
ponent) is reduced by about 10% by the exchange current involving the 7rao(980) exchange 
mechanism. It seems interesting to examine the consequences of these effective coupling 
constants for neutrino-nucleus scatterings. 

To summarize, we wish to reemphasize that neutrino-nucleus reactions play an im- 
portant role in the study of the astrophysical neutrinos. The diversity of possible types of 
reactions presents a great challenge to nuclear physicists. It is hoped that the present re- 
view have succeeded in calling our fellow nuclear physicists' attention to some of the major 
problems facing us in this field. 
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Table 1 Neutrino-deuteron cross sections in units of cm^ 



(MeV) 


vd u'pn 


vd v'pn 


vd ^ e pp 


vd e^nn 


2.0 


O.OOOE+00 


O.OOOE+00 


3.854E-45 


O.OOOE+00 


2.2 


O.OOOE+00 


O.OOOE+00 


8.231E-45 


O.OOOE+00 


2.4 


4.394E-47 


4.362E-47 


1.461E-44 


O.OOOE+00 


2.6 


4.288E-46 


4.253E-46 


2.318E-44 


O.OOOE+00 


2.8 


1.464E-45 


1.451E-45 


3.413E-44 


O.OOOE+00 


3.0 


3.369E-45 


3.334E-45 


4.761E-44 


O.OOOE+00 


3.2 


6.308E-45 


6.236E-45 


6.377E-44 


O.OOOE+00 


3.4 


1.041E-44 


1.028E-44 


8.274E-44 


O.OOOE+00 


3.6 


1.579E-44 


1.557E-44 


1.046E-43 


O.OOOE+00 


3.8 


2.252E-44 


2.219E-44 


1.296E-43 


O.OOOE+00 


4.0 


3.068E-44 


3.021E-44 


1.577E-43 


O.OOOE+00 


4.2 


4.034E-44 


3.967E-44 


1.890E-43 


1.065E-45 


4.4 


5.154E-44 


5.064E-44 


2.237E-43 


4.397E-45 


4.6 


6.434E-44 


6.314E-44 


2.618E-43 


9.832E-45 


4.8 


7.877E-44 


7.722E-44 


3.033E-43 


1.748E-44 


5.0 


9.487E-44 


9.290E-44 


3.485E-43 


2.747E-44 


5.2 


1.127E-43 


1.102E-43 


3.973E-43 


3.993E-44 


5.4 


1.322E-43 


1.292E-43 


4.498E-43 


5.498E-44 


5.6 


1.534E-43 


1.498E-43 


5.060E-43 


7.270E-44 


5.8 


1.765E-43 


1.721E-43 


5.661E-43 


9.319E-44 


6.0 


2.013E-43 


1.961E-43 


6.239E-43 


1.165E-43 


6.2 


2.279E-43 


2.218E-43 


6.904E-43 


1.427E-43 


6.4 


2.563E-43 


2.491E-43 


7.606E-43 


1.719E-43 


6.6 


2.866E-43 


2.782E-43 


8.347E-43 


2.040E-43 


6.8 


3.198E-43 


3.102E-43 


9.127E-43 


2.392E-43 


7.0 


3.541E-43 


3.430E-43 


9.945E-43 


2.774E-43 


7.2 


3.902E-43 


3.776E-43 


1.080E-42 


3.187E-43 


7.4 


4.283E-43 


4.140E-43 


1.170E-42 


3.631E-43 


7.6 


4.683E-43 


4.522E-43 


1.264E-42 


4.106E-43 


7.(S 


5.103E- 13 


1.921E- 13 


1.3G1E-12 


1.G12E-13 
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Table 1 (continued) Neutrino-deuteron cross sections in units of 



(MeV) 


ud — v'pn 


vd — > v'pn 


vd ^ e pp 


— 7 -U 

vd — > e^nn 


8.0 


5.542E-43 


5.339E-43 


1.463E-42 


5.150E-43 


8.2 


6.000E-43 


5.775E-43 


1.569E-42 


5.719E-43 


8.4 


6.479E-43 


6.228E-43 


1.679E-42 


6.320E-43 


8.6 


6.977E-43 


6.700E-43 


1.793E-42 


6.967E-43 


8.8 


7.495E-43 


7.189E-43 


1.911E-42 


7.635E-43 


9.0 


8.033E-43 


7.697E-43 


2.033E-42 


8.336E-43 


9.2 


8.591E-43 


8.223E-43 


2.160E-42 


9.068E-43 


9.4 


9.169E-43 


8.767E-43 


2.290E-42 


9.833E-43 


9.6 


9.767E-43 


9.328E-43 


2.425E-42 


1.063E-42 


9.8 


1.039E-42 


9.908E-43 


2.564E-42 


1.146E-42 


10.0 


1.103E-42 


1.051E-42 


2.708E-42 


1.232E-42 


10.2 


1.168E-42 


1.112E-42 


2.855E-42 


1.322E-42 


10.4 


1.236E-42 


1.176E-42 


3.007E-42 


1.415E-42 


10.6 


1.306E-42 


1.241E-42 


3.163E-42 


1.511E-42 


10.8 


1.378E-42 


1.308E-42 


3.323E-42 


1.610E-42 


11.0 


1.453E-42 


1.377E-42 


3.488E-42 


1.712E-42 


11.2 


1.529E-42 


1.447E-42 


3.657E-42 


1.818E-42 


11.4 


1.607E-42 


1.520E-42 


3.831E-42 


1.927E-42 


11.6 


1.688E-42 


1.594E-42 


4.008E-42 


2.039E-42 


11.8 


1.770E-42 


1.670E-42 


4.191E-42 


2.155E-42 


12.0 


1.854E-42 


1.748E-42 


4.377E-42 


2.273E-42 


12.5 


2.075E-42 


1.950E-42 


4.863E-42 


2.584E-42 


13.0 


2.308E-42 


2.164E-42 


5.378E-42 


2.914E-42 


13.5 


2.555E-42 


2.389E-42 


5.920E-42 


3.265E-42 


14.0 


2.815E-42 


2.625E-42 


6.491E-42 


3.635E-42 


14.5 


3.088E-42 


2.872E-42 


7.090E-42 


4.025E-42 


15.0 


3.375E-42 


3.130E-42 


7.719E-42 


4.435E-42 


15.5 


3.675E-42 


3.399E-42 


8.376E-42 


4.864E-42 


16.0 


3.989E-42 


3.679E-42 


9.063E-42 


5.313E-42 


16.5 


4.316E-42 


3.969E-42 


9.780E-42 


5.780E-42 


17.0 


4.656E-42 


4.271E-42 


1.053E-41 


6.267E-42 


17.5 


5.010E-42 


4.584E-42 


1.130E-41 


6.773E-42 


18.0 


5.378E-42 


4.907E-42 


1.211E-41 


7.297E-42 


18.5 


5.760E-42 


5.241E-42 


1.295E-41 


7.841E-42 


19.0 


6.155E-42 


5.585E-42 


1.381E-41 


8.402E-42 


19.5 


6.564E-42 


5.940E-42 


1.471E-41 


8.983E-42 
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Table 1 (continued) Neutrino-deuteron cross sections in units of 



(MeV) 


vd — > v'pn 


vd — > u'pn 


ud ^ e pp 


vd — > e~^nn 


20.0 


6.987E-42 


6.306E-42 


1.564E-41 


9.581E-42 


20.5 


7.424E-42 


6.682E-42 


1.661E-41 


1.020E-41 


21.0 


7.875E-42 


7.069E-42 


1.760E-41 


1.083E-41 


21.5 


8.339E-42 


7.466E-42 


1.862E-41 


1.149E-41 


22.0 


8.818E-42 


7.873E-42 


1.968E-41 


1.216E-41 


22.5 


9.311E-42 


8.291E-42 


2.077E-41 


1.285E-41 


23.0 


9.818E-42 


8.719E-42 


2.189E-41 


1.355E-41 


23.5 


1.034E-41 


9.157E-42 


2.305E-41 


1.427E-41 


24.0 


1.087E-41 


9.606E-42 


2.423E-41 


1.501E-41 


24.5 


1.142E-41 


1.006E-41 


2.545E-41 


1.577E-41 


25 


1.199E-41 


1.053E-41 


2.671E-41 


1.655E-41 


26 


1.316E-41 


1.150E-41 


2.932E-41 


1.815E-41 


27 


1.439E-41 


1.251E-41 


3.206E-41 


1.981E-41 


28 


1.567E-41 


1.355E-41 


3.494E-41 


2.155E-41 


29 


1.702E-41 


1.464E-41 


3.796E-41 


2.335E-41 


30 


1.842E-41 


1.576E-41 


4.112E-41 


2.521E-41 


31 


1.988E-41 


1.692E-41 


4.443E-41 


2.713E-41 


32 


2.140E-41 


1.812E-41 


4.787E-41 


2.912E-41 


33 


2.298E-41 


1.936E-41 


5.146E-41 


3.117E-41 


34 


2.462E-41 


2.063E-41 


5.520E-41 


3.328E-41 


35 


2.632E-41 


2.194E-41 


5.909E-41 


3.545E-41 


36 


2.808E-41 


2.329E-41 


6.313E-41 


3.768E-41 


37 


2.990E-41 


2.467E-41 


6.731E-41 


3.996E-41 


38 


3.178E-41 


2.609E-41 


7.165E-41 


4.231E-41 


39 


3.372E-41 


2.754E-41 


7.614E-41 


4.470E-41 


40 


3.572E-41 


2.902E-41 


8.079E-41 


4.716E-41 


41 


3.778E-41 


3.055E-41 


8.559E-41 


4.966E-41 


42 


3.991E-41 


3.210E-41 


9.055E-41 


5.222E-41 


43 


4.209E-41 


3.368E-41 


9.567E-41 


5.484E-41 


44 


4.434E-41 


3.530E-41 


1.009E-40 


5.750E-41 


45 


4.665E-41 


3.695E-41 


1.064E-40 


6.021E-41 


46 


4.901E-41 


3.864E-41 


1.120E-40 


6.298E-41 


47 


5.144E-41 


4.035E-41 


1.177E-40 


6.579E-41 


48 


5.393E-41 


4.210E-41 


1.237E-40 


6.865E-41 


49 


5.649E-41 


4.387E-41 


1.297E-40 


7.156E-41 
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Table 1 (continued) Neutrino-deuteron cross sections in units of 



E„ (MeV) 


ud ^ V pn 


— 7 — / 

vd ^ V pn 


vd ^ e pp 


— 7 4- 

vd e^nn 


50 


5.910E-41 


4.568E-41 


1.360E-40 


7.451E-41 


51 


6.177E-41 


4.751E-41 


1.424E-40 


7.751E-41 


52 


6.451E-41 


4.937E-41 


1.490E-40 


8.055E-41 


53 


6.731E-41 


5.126E-41 


1.557E-40 


8.364E-41 


54 


7.017E-41 


5.318E-41 


1.626E-40 


8.676E-41 


55 


7.309E-41 


5.513E-41 


1.697E-40 


8.993E-41 


60 


8.8dOE-41 


6.525E-41 


2.076E-40 


1.064E-40 


65 


1.056E-40 


7.599E-41 


2.498E-40 


1.237E-40 


70 


1.241E-40 


8.728E-41 


2.961E-40 


1.419E-40 


75 


1.440E-40 


9.906E-41 


3.467E-40 


1.607E-40 


80 


1.654E-40 


1.113E-40 


4.014E-40 


1.801E-40 


85 


1.881E-40 


1.239E-40 


4.603E-40 


2.001E-40 


90 


2.121E-40 


1.368E-40 


5.232E-40 


2.205E-40 


95 


2.373E-40 


1.500E-40 


5.900E-40 


2.412E-40 


100 


2.637E-40 


1.634E-40 


6.606E-40 


2.623E-40 


110 


3.199E-40 


1.908E-40 


8.127E-40 


3.049E-40 


120 


3.801E-40 


2.186E-40 


9.784E-40 


3.481E-40 


130 


4.437E-40 


2.465E-40 


1.156E-39 


3.914E-40 


140 


5.103E-40 


2.744E-40 


1.345E-39 


4.346E-40 


150 


5.793E-40 


3.020E-40 


1.543E-39 


4.776E-40 


160 


6.501E-40 


3.292E-40 


1.750E-39 


5.201E-40 


170 


7.224E-40 


3.559E-40 


1.963E-39 


5.623E-40 
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